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ABSTRACT 

In Paper III of our series “A Uniform Analysis of the Ly-ct forest at z — 0 — 5”, 
we presented a set of 270 quasar spectra from the archives of the Faint Object 
Spectrograph (FOS) on the Hubble Space Telescope (HST). A total of 151 of 
these spectra, yielding 906 lines, are suitable for using the proximity effect signa¬ 
ture to measure J(v o), the mean intensity of the hydrogen-ionizing background 
radiation held, at low redshift. Using a maximum likelihood technique and the 
best estimates possible for each QSO’s Lyman limit flux and systemic redshift, 
we find J(u Q ) = 7.6l|;o x 10” 23 ergs s” 1 cm -2 Hz” 1 sr” 1 at 0.03 < z < 1.67. This 
is in good agreement with the mean intensity expected from models of the back¬ 
ground which incorporate only the known quasar population. When the sample 
is divided into two subsamples, consisting of lines with z < 1 and z > 1, the val¬ 
ues of J(v 0 ) found are 6.5+ 38 6 x 10” 23 ergs s” 1 cm” 2 Hz” 1 sr” 1 , and I.CUq; 8 x 10” 22 
ergs s” 1 cm” 2 Hz” 1 sr” 1 , respectively, indicating that the mean intensity of the 
background is evolving over the redshift range of this data set. Relaxing the 
assumption that the spectral shapes of the sample spectra and the background 
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are identical, the best fit HI photoionization rates are found to be 6.7 x ICC 13 
s^ 1 for all redshifts, and 1.9 x 10~ 13 s” 1 and 1.3 x 10~ 12 s _1 for z < 1 and 
z > 1, respectively. The inclusion of blazars, associated absorbers, or damped 
Ly-ct absorbers, or the consideration of a A CDM cosmology rather than one in 
which 11 a = 0 has no significant effect on the results. The result obtained using 
radio loud objects is not significantly different from that found using radio quiet 
objects. Allowing for a variable equivalent width threshold gives a consistently 
larger value of J(v o) than the constant threshold treatment, though this is found 
to be sensitive to the inclusion of a small number of weak lines near the QSO 
emission redshifts. This work confirms that the evolution of the number density 
of Ly-Q: lines is driven by a decrease in the ionizing background from z ~ 2 to 
z ~ 0 as well as by the formation of structure in the intergalactic medium. 

Subject headings: diffuse radiation— intergalactic medium—quasars: absorption 
lines 


1. Introduction 

The spectra of quasars show a “forest” of absorption lines blueward of the Ly-o emission 
line (Lynds 1971, Sargent et al. 1980, Weymann, Carswell, & Smith 1981). Observational 
and theoretical work in recent years has shown that most of this absorption can be attributed 
to neutral hydrogen in galaxies and large-scale structure along the line of sight (Cen et al. 
1994, Lanzetta et al. 1995,1996, Zhang et al. 1995, Hernquist et al. 1996, Miralda-Escude 
et al. 1996, Bi & Davidsen 1997, Chen et al. 1998, Theuns et al. 1998, Ortiz-Gil et al. 
1999, Impey, Petry, & Flint 1999, Dave et al. 1999, Bryan et al. 1999). In aggregate, QSO 
spectra show an increasing line density with increasing redshift such that dN/dz oc (1. + z) 1 
(Sargent et al. 1980, Weymann, Carswell, & Smith 1981, Young et al. 1982, Murdoch et al. 
1986, Lu, Wolfe, & Turnshek, 1991, Bechtold 1994, Kim et al. 1997). But the line density 
within an individual quasar spectrum decreases with proximity to the Lj-a emission line 
(Weymann, Carswell, & Smith 1981, Murdoch et al. 1986). This is generally thought to be 
due to enhanced ionization of neutral hydrogen in the vicinity of the quasar due to ionizing 
photons from the quasar itself. This “proximity effect” can be used to measure the mean 
intensity of the UV background, denoted J{y 0 ) (Carswell et al. 1987, Bajtlik, Duncan, & 
Ostriker 1988, hereafter BDO). 

J{y o) has been measured at z > 1.7 by a variety of authors (BDO, Lu, Wolfe, & 
Turnshek 1991, Giallongo et al. 1993,1996, Bechtold 1994, Williger et al. 1994, Cristiani 
et al. 1995, Fernandez-Soto et al. 1995, Lu et al. 1996, Savaglio et al. 1997, Cooke et al. 
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1997, Scott et al. 2000b). The results, summarized in Paper II of this series (Scott et al. 
2000b), are in general agreement with the predictions of models of the UV background which 
integrate the contribution from known population of quasars and include reprocessing effects 
in an inhomogeneous intergalactic medium (Haardt & Madau 1996, hereafter HM96, Fardal 
et al. 1998). In Paper II, the mean intensity of the ionizing background was found to be 
7.01^4 x 10' 22 ergs s' 1 cm -2 Hz -1 sr' 1 at z ~ 3. The decline of the quasar space density 
from z ~ 2 to the present is expected to drive a corresponding decline in the intensity of the 
UV background. Kulkarni & Fall (1993, hereafter KF93) measured J{v 0 ) ~ 6 x 10 -24 ergs 
s' 1 cm -2 Hz' 1 sr' 1 at z ~ 0.5 from a subset of the now complete HST Quasar Absorption 
Line Key Project sample presented by Bahcall et al. (1993). 

Much of this previous work has relied upon the technique for measuring J(z/ 0 ) outlined 
by BDO. This technique requires the entire sample of absorption lines to be binned according 
to the ratio of the quasar flux at the physical position of the absorber to the background 
flux. This is done for several initial guesses of the background intensity; and the value that 
gives the lowest y 2 between the binned data and the ionization model is chosen as the best 
fit J(u 0 ). However, this is not the optimal technique to use at low redshift where absorption 
line densities are low. KF93 developed a maximum likelihood technique to address this issue 
and used it in their measurement of J(z/ 0 ) at z ~ 0.5. However, their measurement was based 
upon a sample of only 13 QSOs and less than 100 lines, and has correspondingly large error 
bars. In addition, the value these authors find is lower than the predictions of the models 
of Haardt & Madau (1996), though consistent within the uncertainties, as shown in Figure 
13 of Paper II and in Figure 11 of this paper. Given the importance of the value of the HI 
ionization rate to the hydrodynamical evolution of the low redshift universe, performing this 
measurement with a much larger line sample is worthwhile. 

The low redshift hydrodynamic simulations of Theuns et al. (1998) and Dave et al. 
(1999) indicate that the evolution of the ionizing background is the primary driver behind 
the change of character of the Ly-a forest from high redshift to low redshift, specifically, the 
break in the number distribution of Ly-a lines at z — 1.7 (Morris et al. 1991, Bahcall et al. 
1991, Weymann et al. 1998). The growth of structure pulling gas from low density regions 
into high density regions also contributes to this and other attributes of the evolution of the 
Ly-a forest. 

Shull et al. (1999) estimate the local ionizing background including contributions to the 
background from starburst galaxies as well as Seyferts and QSOs. Their models include a 
treatment of the opacity of the low redshift Ly-a; forest using information drawn from recent 
observational work (Weymann et al. 1998, Penton et al. 2000b). They find that starbursts 
and AGN could contribute approximately equally to the ionizing background at low redshift, 
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each ~ 1.0 x 10” 23 ergs s” 1 cm” 2 Hz -1 sr” 1 . 

The full HST/FOS archival data set is presented in Paper III and can also be found at 

http: //lithops . as. arizona.edu/~jill/QuasarSpectra or 

http://hea-www.harvard.edu/QEDT/QuasarSpectra, 

so we describe the data used in this paper only briefly in §2. In §3 and §4, we discuss our 
treatment of two parameters of each sample object which are integral to the proximity effect 
analysis, systemic redshifts and Lyman limit fluxes. We outline the proximity effect analysis 
in §5 and we present onr results in §7. We discuss the recovery of J{v o) from simulated Ly-ct 
forest spectra in §6.1. A comparison of the results from radio loud and radio quiet QSOs 
is given in §7.1. The maximum likelihood solutions for J(z/ 0 ) found when allowing for an 
equivalent width threshold that varies across each sample spectrum are discussed in §6.3. 
Solutions for the HI ionization rate are given in §6.2. The effect of a non-zero 12 a on our 
calculations is discussed in §7.2, and the effect of the ionizing background on the Ly-a forest 
line density is discussed in §7.3. We provide comparisons with previous observational work 
on the low redshift UV background in §7.4 and with models of this background in §7.5. A 
discussion of possible systematic effects on this analysis is given in §7.6. We conclude with 
a summary of the results in §8. 


2. Data Sample 

The reduction of the FOS data is described in Paper III. Table 1 lists the objects used in 
the proximity effect analysis along with the object’s redshift and classification in the NASA 
Extragalactic Database. 

For the reasons outlined in Scott et al. (2000b) we have removed from the full FOS 
sample of Paper III the spectra of quasars known to be lensed, as well as those that show 
damped Ly-a absorption, associated absorption, or broad intrinsic absorption. For onr pri¬ 
mary proximity effect sample, we also remove objects classified as blazars (BL Lacs and 
optically violent variables) on the grounds that their continua are highly variable. How¬ 
ever, we also perform the proximity effect analysis with associated absorbers, damped Ly-a 
absorbers, and blazars included in order to determine if they affect the results obtained. 

As discussed in Paper III, objects observed only in the period before the COSTAR 
upgrade to the HST optics and with the A-l FOS aperture are particularly subject to 
irregular line spread functions. We have omitted those data from this analysis as well. The 
distributions in redshift of the QSOs and absorption lines used in this paper are shown in 
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Figure 1. 


3. Systemic Redshifts 

QSO redshifts based on the Ly-a emission line have been shown to be blueshifted from 
the systemic redshift based on narrow emission lines by up to ~ 200 km s -1 . Generally, the 
forbidden OII1 doublet at 4959, 5007 A is taken to be the most reliable indicator of the QSO 
systemic redshift; though other lines such as Mg If AA2796,2803 and Balmer lines have been 
shown to trace the systemic redshift as well, with some spread. (Zheng & Sulentic 1990, 
Tytler & Fan 1992, Laor et al. 1994,1995, Corbin & Boroson 1996) However, the results of 
Nishihara et ah (1997), M c Intosh et ah (1999), and Scott et ah (2000b) indicate that in fact 
H (3 may not reflect the systemic redshift of high redshift QSOs. 


3.1. Observations 

Spectra of the emission lines H/J, [OIII] A5007, or Mg 11 were obtained for several objects 
in our total proximity effect sample. The observations were carried out on the nights of 19 
December 1995, 14 January 1996, 20 and 21 April 1996, 12 and 13 December 1996, and 2 
February 1997. These observations are summarized in Table 2. 

The 19 December 1995 and 13 December 1996 observations were made using the 1.5 
meter Tillinghast telescope at the Fred Lawrence Whipple Observatory using the FAST 
spectrograph (Fabricant et al. 1998) and a thinned Loral 512x2688 CCD chip (gain = 1.06, 
read noise = 7.9 e“) binned by a factor of 4 in the cross-dispersion direction. Observations 
were made using a 300 lines mm -1 grating blazed at 4750 A and a 3" slit. These spectra 
cover a wavelength range of 3660-7540 A. This is listed as set-up (1) in Table 2. 

The January, April, and December 10 and 12, 1996 observations were made using the 
Steward Observatory Bok 90 inch telescope using the Boiler and Chivens Spectrograph with 
a 600 1 mm -1 grating blazed at 6681 A in the first order, a 1.5" slit, and a 1200 x 800 
CCD array with a gain of 2.2 e~ ADU” 1 and a read noise of 7.7 e _ , binned lxl. For 
the January 1996 observations, the data were obtained with one of two grating tilts, one 
resulting in wavelength coverages of 3600-5825 A and 6870-9140 A. For the April 1996 data, 
the wavelength ranges were 4140-6370 A and 5280-7550 A. Two grating tilts were also used 
for the December 1996 data, giving wavelength coverages of 4500-6700 A and 5610-7860 A. 

The spectrum of one object, 0827+2421, was obtained on 15 February 1997 at the 
Multiple Mirror Telescope with the Blue Channel Spectrograph, a 2 " slit, the 3K x IK CCD 
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array, and the 800 1 mm 1 grating blazed at 4050 A with spectral coverage of 4365-6665 A. 

The spectra are shown in Figure 2 and the lines used for redshift measurements are 
labeled. 


3.2. Measurements 

Taking a simple cursor measurement of each line centroid, we find a mean [OIII]-Balmer 
line Av of -30 ± 1010 km s -1 for 31 objects and a mean [0I11]-Mg II Av of 58 ± 576 km s -1 
for 31 objects. The mean blueshift of the Ly-a emission line with respect to [OIII] is 289 
± 727 km s -1 based on 51 measurements. The redshifts measured for each object in our 
sample are shown in Table 1; and the results are shown in Figure 3. Gaussian fits to the 
lines give similar results. 

We therefore treat both Balmer lines and Mg II in addition to [OIII] as good systemic 
redshift indicators for these low redshift objects. In the case of a QSO for which we have 
only a Ly-a emission line measurement of the redshift, we add 300 km s -1 to this value to 
estimate its systemic redshift. 


4. Lyman Limit Fluxes 

Our method for estimating Lyman limit fluxes for each QSO is the same as that described 
in Paper II. For objects with spectral coverage between the Ly -a and CIV emission lines, 
we extrapolate the flux from 1450 A in the quasar’s rest frame to 912 A using /„ ~ u~ a 
and a spectral index a measured primarily from the spectral region between the Ly-a and 
C IV emission lines. Figure 4 shows the FOS spectra for which these fits were made along 
with the power law fits themselves. In some cases, a is poorly constrained from these fits, 
especially if there was little spectral coverage redward of Ly-a emission in the data. If 
another measurement of the spectral index was available in the literature for these objects, 
we used it; otherwise, we used our measurement. 

Table 3 lists the Lyman limit flux for each object in this proximity effect sample and 
either a) the flux at 1450 A, or some other appropriate wavelength free of emission features, 
measured from the FOS data, or b) a directly measured Lyman limit flux and the reference. 
If available from the extracted archive data, red spectra and the fits to them are presented for 
objects which were observed only with pre-COSTAR FOS and A-l aperture, though these 
data were not subsequently used for any Ly-a forest studies. See Table 4 of Paper III. 
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In Figure 5, we show QSO Lyman limit luminosities versus emission redshift for this 
HST/FOS sample combined with the high redshift objects presented in Papers I and II. Only 
at the lowest redshifts is there any trend of luminosity with redshift. 


5. Analysis 


The distribution of Ly-a lines in redshift and equivalent width is given by: 


d 2 M 

dzdW 


Ap 

W* 


(1 + z) 1 exp 



( 1 ) 


The distribution in redshift and HI column density, N, is: 

r) 2 Kf 

— = AN -»(i + z y. ( 2 ) 

The parameter 7 is the redshift distribution parameter. The quantities W* in Equ. 1 and 
f3 in Equ. 2 are the line rest equivalenth width and column density distribution parameters, 
respectively. The quantities Ap and A are normalizations. 


The BDO method for measuring J(vp) consists of binning all lines in the sample in the 
parameter u(z), the ratio of QSO to background Lyman limit flux density at the physical 
location of the absorber: F Q (v 0 ) / (An J (up)) for various values of J(vp ). The value of J(vp) 
that results in the lowest y 2 between the binned data and the ionization model, 

'A- = _4 0 (1 + *) 7 [1 +u<z)]- ( ' , ~ 1) , (3) 

is considered to be the optimal value. This ionization model follows from the assumption 
that the column densities of lines are modified by the presence of the QSO according to 


N oc N 0 (l + u(z)) \ 


(4) 


where N 0 is the column density a given line would have in the absence of the QSO. The la 
errors are found from Ay 2 = 8.18 for 7 degrees of freedom (Press et al. 1992). 


The value of uj(z) for each line in a given sample depends not only upon the value of 
J(vp) assumed, but also on the cosmological model, as 


F Q (v 0 ) 


L(vp) 

Awrl(z) 


(5) 


/M 

(1 + Am) ’ 


and 


L(v 0 ) = 4jrt/|0) 


( 6 ) 



where tl{z) is the luminosity distance of an individual absorber from the QSO and d^z) 
luminosity distance to the QSO from the observer. The luminosity distance between two 
objects at different redshifts can be calculated analytically for cosmological models in which 
Oa = 0. We return to this point in Section 7.2 below. 


If the proximity effect is indeed caused by enhanced ionization of the IGM in the vicinity 
of QSOs, one may expect to observe a larger deficit of lines relative to the Ly-a forest near 
high luminosity QSOs than near low luminosity QSOs. In Figure 6 (a), we plot the fractional 
deficit of lines with respect to the number predicted by Equ. 1 versus distance from the QSO 
for this HST/FOS sample combined with the high redshift objects observed with the Multiple 
Mirror Telescope (MMT) presented in Papers I and II. We divide our QSO sample into high 
and low luminosity objects at the median Lyman limit luminosity of the combined MMT 
and HST/FOS sample, log(L 912 ^) ~ 31. High luminosity objects show a marginally more 
pronounced proximity effect than low luminosity objects: 4.9a for QSOs with log(L 912 > 31 
versus 3.2 a for QSOs with log(L 912 ^) < 31. In panel (b), we plot the line deficit within 2 
hyd Mpc as a function of log(L 912 ^). The lack of a significant difference in the line deficit 
between high and low luminosity QSOs may indicate the presence of clustering, if absorption 
features cluster more strongly around more luminous QSOs with deeper potential wells. We 
will address the issue of clustering further below. 


The BDO method of measuring the background can result in poor statistics at low 
redshift due to the low line density in the low redshift Ly-a forest. We will quote results 
from this method, but we will generally the maximum likelihood method for measuring J{v 0 ) 
as presented by KF93, which consists of constructing a likelihood function of the form 


L 


n/(^)n-p[- 

a q 



' N q . 

min 


f(N,z)dN], 


(7) 


where 

/(IV, z ) = AN~ P { 1 + Z y\l + ( 8 ) 

and the indicies a and q denote sample absorption lines and quasars, respectively. Using the 
values of 7 and A 0 from a separate maximum likelihood analysis on the Ly-a forest excluding 
regions of the spectra affected by the proximity effect (Dobrzycki et al. 2001, Paper IV), 
and a value of (3 from studies with high resolution data, eg. /3 = 1.46 from Hu et al. (1995), 
the search for the best-fit value of J{y 0 ) consists of finding the value that maximizes this 
function, fixing the other parameters. 


If the line density is low throughout a single Ly-a: forest spectrum, it becomes difficult to 
distinguish any proximity effect, even in a large sample of spectra. The absence of a proximity 
effect in this model formally translates into the limit J{v 0 ) —> 00 because in this scenario, 
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the QSO has no additional effect on its surroundings and therefore generates no relative 
line underdensity. The errors quoted in the values of log[J(z/ 0 )] are found from the fact 
that in solving for log[J(z/ 0 )] alone, the logarithm of the likelihood function, — 21 n(L/L max ), 
is distributed as y 2 with one degree of freedom. In the case of an ill-defined solution, the 
likelihood function is very broad and the formal error approaches infinity. If a proximity effect 
is weak but not absent in the data, a maximum likelihood solution is sometimes possible, 
but with no well-defined 1 a upper limit on log[ J(z/ 0 )]- In other words, if an upper limit of 
infinity is quoted, the data cannot rule out the nonexistence of a proximity effect to within 
la confidence. 

Using a constant equivalent width threshold results in the loss of a large amount of 
spectral information. In the case of a large equivalent width threshold, of course, many weak 
lines are discarded; and in the case of a small threshold, regions of spectra where the signal- 
to-noise ratio (S/N) does not permit the detection of lines all the way down to the specified 
threshold are lost and only the highest S/N spectral regions are used. The technique of 
measuring the statistics 7 and W* has been expanded to allow for a threshold that varies 
with S/N across each QSO spectrum (Bahcall et al. 1993,1996, Weymann et al. 1998, Scott 
et al. 2000a). We will use this variable threshold information to measure J{v 0 ) as well. 


6. Results 


The results of this analysis are given in Table 4. 

Before we begin the discussion of the results, some words about the normalization values 
listed in Table 4 are in order. In the BDO method for measuring J{v 0 ), lines are binned 
in u(z) and compared to the ionization model given by Equ. 3, for an assumed value of (3. 
In this case, the normalization listed in Table 4 is the parameter in Equ. 1, found from the 
number of lines in the sample and the maximum likelihood value of 7 : 


A 0 = A 0 exp 


Wi im \ 

it; ) 



z q 

^min 


dz (1 + z)' 


-1 


5 


(9) 


where A f is the total number of lines observed with rest equivalent width greater than Wu m , 
the limiting equivalent width of the line sample. For the maximum likelihood solutions for 
J{y 0 ), we convert line equivalent widths to column densities using the Ly-a curve of growth 
and an assumed value of b , the characteristic Doppler parameter of the lines. As we will 
demonstrate, different values of f3 and b have only a small effect on the value of J{y 0 ) found. 
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The normalization is given by 




dNN-^l+z) 



-1 


5 


( 10 ) 


where N^- n (z) is the limiting column density across each QSO spectrum corresponding to a 
limiting equivalent width. This quantity can be held constant, as in the BDO method, or 
it can be allowed to vary across each QSO spectrum. In both of these formulations for the 
normalization, a proximity region around the QSO is neglected and that proximity region is 
either defined by a velocity cut, eg. z em - 3000 km s^ 1 , or by a cut in uj(z), eg. c u(z) = 0.2. 

We also use the standard BDO method to find log[ J(z/ 0 )] = — 22 . 041^11 and —22.06^62 
for equivalent width thresholds of 0.32 and 0.24 A respectively. Figures 7(a) and (d) illustrate 
the x 2 of the binned data compared to the BDO ionization model as a function of assumed 
J{y 0 ) for these two thresholds. The BDO ionization model is expressed in terms of the 
number of lines per coevolving coordinate: 


dN 

dXj 


Ao(l + a;) 1 \ 


( 11 ) 


where X 7 = f (1 + z) 7 dz. This y 2 curve is very broad, which is reflected in the large error 
bars and indicates the difficulty in isolating the optimal mean intensity of a weak background 
using this technique. Figures 8 (a) and (d) show the binned data and the ionization model 
for the values of J(u 0 ) listed above, those that give the lowest y 2 between the binned data 
and the model, ie. the minima of the curves in Figures 7(a) and (d). 

We executed the maximum likelihood search for J{v 0 ), using two different fixed equiva¬ 
lent width thresholds, 0.24 A and 0.32 A as well as for the case of a variable threshold across 
all the spectra. The uncertainty in 7 does not translate directly into a large uncertainty in 
J{v 0 ). Changing the value of 7 alters the maximum likelihood normalization, A, according 
to Equ. 10 . From the sample of lines with rest equivalent widths greater than 0.32 A we find 
log [J(u 0 )} = -22.11±%% for 7 = 0.82 ± 0.29. Varying 7 by ±la gives log[J(i/ 0 )] = - 22.21 
and — 22.00 with similar uncertainties. 


The data used here are not of sufficient resolution to fit Voigt profiles to the absorption 
features and derive HI column densities and Dopper parameters. We therefore choose to fix 
the values of /3 and b to those found from work on high resolution data, rather than allow 
them to freely vary in our analysis. For the 0.32 A fixed equivalent width threshold, we 
tested several pairs of values of (/ 3 , b) where b is in km s -1 : (1.46,35) and (1.46,25) where the 
value of (3 is taken from Hu et al. (1995); as well as (1.45,25) and (1.70,30) found from low 
redshift Ly-cx forest spectra taken with the Goddard High Resolution Spectrograph (GHRS) 
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on HST by Penton et al. (2000a,b). In addition, Dave & Tripp (2001) have found some 
evidence for (3 increasing to 2.04 at z < 0.3 from high resolution echelle data from the Space 
Telescope Imaging Spectrograph aboard the HST. We test this value as well. The likelihood 
functions for the maximum likelihood solutions listed in rows 2-6, 8-12, 14, and 18 of Table 4 
are shown in Figure 9. The binned data and ionization models are plotted in Figure 10. 
The values of J{v o) derived for these various pairs of values of [3 and b are not significantly 
different from one another, though the results in Table 4 indicate that varying f3 has a larger 
impact on the inferred J(u 0 ) than does varying b. The solution for (3 = 2.04 differs from the 
(3 = 1.46 solution by ~ la. In the analysis that follows, we adopt the values 1.46 and 35 
km s -1 . 

The models of Haardt & Madau (1996) predict that the UV background arising from 
QSOs drops by over an order of magnitude from z = 2.5 to z = 0. We therefore divide 
the sample into low and high redshift subsamples at z — 1 and use both the BDO method 
and the maximum likelihood method for finding J{v o). These results, also listed in Table 4, 
confirm some evolution in J{v o), though not at a high level of significance. For the BDO 
solutions, we find log[./(z/ 0 )] at z < 1 is equal to — 22 . 87 to;g 2 and log[J(z/ 0 )] at z > 1 is 
equal to — 22.02l ( j l '33 5 - The restrictive la upper limit for log[J(^o)] at z > 1 arises from the 
steeply rising y 2 as a function of log[,7(// 0 )] shown in Figure 7. This, in turn arises from 
the single line in the highest log(u;) bin moving to the next bin for larger values of J{v o), 
resulting in a drastic change in the y 2 with respect to the photoionization model. We do not 
consider this to be a reliable indicator of the uncertainty in J{v 0 ) at z > 1. The maximum 
likelihood technique gives more robust estimates of the uncertainties. From this analysis, we 
find log[J(z/ 0 )] at z < 1 is found to be — 22 .181^61, while at z > 1 it is -21.98lo;^. These 
results are shown in Figures 11(a) and 16. 

Including associated absorbers, clamped Ly-a absorbers, or blazars in the proximity 
effect analysis appears to have little effect on the results. One might expect associated 
absorbers to reduce the magnitude of the observed proximity effect and hence cause J(uq) 
to be overestimated. The value found including the 45 associated absorbers in our sample is 
indeed larger, log[J(z/ 0 )]= — 21 . 74 ^ 39 , versus log[J(z/ 0 )]= — 22 . 11 ^ 40 , but not significantly 
so. Likewise, if the intervening dust extinction in damped Ly-a absorbers is significant, 
including these objects in our analysis could cause us to overestimate the magnitude of the 
proximity effect and hence underestimate J(u o). However, the inclusion of these 7 objects 
only negligibly reduces the value of J{v 0 ) derived. QSO variability on timescales less than 
~ 10 5 years would be expected to smooth out the proximity effect distribution (BDO). 
However, the inclusion of 6 blazars in the sample, all at z < 1, resulted in no discernible 
change in J{y o). The sample used in the analysis of HI ionization rates discussed below 
includes all of these objects. 
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For each solution, we calculate the y 2 with respect to the ionization model expressed by 
Equ. 3, and the probability that the observed y 2 will exceed the value listed by chance for a 
correct model, Q x 2 (Press et al. 1992). We also execute a Kolmogorov-Smirnov (KS) test for 
each solution. The KS test provides a measure of how well the assumed parent distribution 
of lines with respect to redshift, given by Equ. 8 , reflects the true redshift distribution of 
lines (cf. Murdoch et ah 1986, Press et al. 1992). The KS probability, Qks, indicates the 
probability that a value of the KS statistic larger than the one calculated could have occurred 
by chance if the assumed parent is correct. The KS probability associated with each solution 
for J(y o) is listed in column 10 of Table 4. 


6.1. Simulations 

We tested our maximum likelihood methods, including our treatment of the variable 
equivalent width thresholds by running our analysis on a simulated data set. Each of the 151 
spectra in this simulated data set had a redshift equal to that of an object in our data set. All 
objects including those showing associated absorption, damped Ly -a absorption, or blazar 
activity are included in this simulated set. Each spectrum is created using a Monte Carlo 
technique by which lines are placed in redshift and column density space according to Equ. 2. 
A background of known mean intensity modifies the column densities of the lines according to 
the BDO formulation given by Equ. 4. The same analysis done on the data, consisting of the 
line-finding algorithm and the maximum likelihood searches for 7 and J{v 0 ), is then used on 
the simulated spectra in order to recover the input J[y 0 ). Three different values of log[ J(z/ 0 )] 
are input, -21, -22, and -23, and the results are listed in Table 5. In order to understand the 
possible range of recovered log[J(z/ 0 )], we repeated the input log[J(z/ 0 )]= —22 simulation in 
the constant threshold case nine additional times, resulting in J[v 0 ) = 2.91 ± 1.67 x 10 ~ 22 
ergs s _1 cm -2 Hz -1 sr -1 . In addition, since we observe the background to evolve with redshift 
from z = 1 to z = 0, we implement a model in which J{y 0 ) varies as a power law in (1 + z) 
over the redshift range of the data. This relationship is defined by the best fit to a power 
law variation of J{y 0 ) with redshift: log[J(z/ 0 )] = 0.017 log(l + z) — 21.87. We recover this 
using both the constant threshold and the variable threshold analyses, at all redshifts and 
at z < 1 and z > 1 separately. The results of this exercise are shown in Table 5 and in 
Figure 12 . 

These simulation results indicate that both the constant and variable threshold analyses 
can overestimate the background by up to a factor of 3-5, though the uncertainties for the 
variable threshold solutions are consistently lower, as a factor of ~ 2 more lines are used in 
these solutions. We separated the first of the input log[J(z/ 0 )]= —22 simulated data samples 
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into high and low redshift subsamples at z = 1, in order to determine if the change in J(v o) 
as a function of redshift could be falsely introduced in a case there the input background is 
constant with redshift. For both the constant and variable threshold treatments, this is not 
the case. The value found for the low redshift subsample is actually larger than the value 
found for the high redshift subsample in both treatments. 

In the case of the varying input log[J(z/ 0 )], the values recovered for the high redshift 
subsample and for the entire redshift range of the data are overestimates. The slope of 
the linear relationship between log[ J(z/ 0 )] and log(l + z) is quite small, 0.017, resulting in 
a variable input log[./(Vo)] that is actually nearly constant with redshift. The solution for 
z < 1 matches the input well for both the constant and variable threshold cases. At z > 1, 
the variable threshold solution overestimates the input by a larger factor, ~3, or 1.6er, than 
does the constant threshold solution, ~2, or less than la. 

In Paper II, we argued that curve-of-growth effects are likely to come into play in the 
proximity effect analysis and to play a larger role for cases in which ./(Vo) is large and 
the proximity effect signature is small. Here we find that the input J(v o) is recovered 
most effectively by the constant and variable threshold cases for the largest input value of 
log[J(z/ 0 )], —21. However, nearly every case tested with these simulations results in a value 
of J{y o) larger than the input value, especially when a variable equivalent width threshold is 
used. The only case where the difference is significant is the input log[J(Vo)] = —23, variable 
threshold case. The recovered log[J(i/o)], -22.47, is 4 a larger than the input. We will return 
to the discussion of the variable threshold in Section 6.3 below. 


6.2. HI Ionization Rate 


As described in Paper II, solving for the HI ionization rate, 


r 



AT[J{y)a HI {y) 

--- au s 

hv 




( 12 ) 


instead of J{v 0 ) avoids the assumption that the spectral indicies of the QSOs and the back¬ 
ground are identical. We modified our maximum likelihood code to use the values of a for 
each QSO listed in Table 3 to measure this quantity and the results are listed in Table 6. 
For objects with no available measured value of a, we use a = 2.02, the extreme ultraviolet 
spectral index measured from a composite spectrum of 101 HST/FOS QSO spectra by Zheng 
et al. (1997). The result for lines above a constant 0.32 A rest equivalent width threshold 
is log(r) = — 12.17tg;4Q. This result is not substantially changed if we instead use a = 1.76, 
the value found from a composite of 184 QSO spectra from HST/FOS, GHRS, and STIS 
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by Tclfer et al. (2001), giving log(r) = —12.25^35• We also find little change in the result 
if we assume a = 2.02 or a = 1.76 for all QSOs. The variable threshold data result in a 
high HI ionization rate, and this is discussed further in the following section. The constant 
threshold result is plotted in Figure 13. Evolution in the UV background is more apparent 
in the HI ionization rate than in the solutions for J(y o). The result at z > 1 is 6.5 times 
larger than that at z < 1. The values of J{v o) implied by these solutions for T and a global 
source spectral index a s = 1.8 are also listed in Table 6. 

We also parametrize the evolution of the HI ionization rate as a power law: 

r(z) = A pl (i + z) B * (13) 


and solve for the parameters A p i and B p \ in both the constant and variable threshold cases. 
The values we find are shown as the dashed line in Figure 13 also listed in Table 6. 

HM96 parametrize their models of the HI ionization rate with the function: 


r(z) = H hm (1 + -) Shm 


exp 


-\z - z c 

S 


(14) 


We combine our data set with that of Scott et al. 2000b to solve for the parameters A H m, 
Bum, z c , and S. We find (A H m, -E>hm, z c , S)= (7.6 x 1CT 13 , 0.35, 2.07, 1.77) for /?=1.46 and 
(Hhm, Bum, z c , S) = (3.2 x 10 -13 , 1.45, 2.13, 1.42) for (3=1.7, while the parameters found by 
HM96 for q 0 = 0.5 are (6.7 x 1CT 13 , 0.43, 2.30, 1.95). These results are also represented by 
the solid curves in Figure 13, while the HM96 parametrization is shown by the dotted line 
for comparison. 


6.3. Variable Equivalent Width Threshold 

The variable threshold analysis yielded some unexpected results. As seen in the majority 
of the simulations, the values of J{y o) found were consistently larger than the values found 
using a constant equivalent width threshold, indicating that the inclusion of weaker lines 
suppresses the proximity effect. This is to be expected if clustering is occurring (Loeb & 
Eisenstein 1995), which in itself is to be expected to be more prominent at low redshift than 
at high redshift. However, the suppression of the proximity effect by the inclusion of weak 
lines is somewhat counterintuitive from the perspective of the curve of growth. Most of the 
lines included in a constant threshold solution are on the flat part of the curve of growth. 
Therefore, though the ionizing influence of the quasar may be translated directly into a 
change in the HI column density, as predicted by the BDO photoionization model, this will 
not necessarily result in a corresponding change in the line equivalent width. The solution 
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for z < 1 is nearly a factor of 3 larger than the the solution found in the case of a constant, 
0.32 A equivalent width threshold. The solution for z > 1 is a factor of ~ 6 larger than the 
constant threshold solution, with no well-defined la upper limit due to the flattening of the 
likelihood function towards high J(u 0 ) This likelihood function for the total sample shows 
two peaks, the most prominent at log[J(z/ 0 )]= —20.82, the solution listed in Table 4, and a 
secondary peak at log[J(z/ 0 )] ~ —18.4. 

This behavior is also exhibited, even more dramatically, in the solutions for the HI 
ionization rate, as discussed above. We conducted a jackknife resampling experiment (Babu 
& Feigelson 1996, Efron 1982) to determine the source of these likelihood function peaks at 
large log(T), or log[ J(u 0 )\. 

Two objects, 0743-6719 (z em = 1.508) and 0302-2223 (z em = 1.402), are found from 
jackknife experiments to produce all of this effect. In the jackknife experiment, we perform 
the maximum likelihood calculation of J{v o) N times, where N is the number of objects 
in the high redshift subsample. In each calculation, one object from the total sample is 
removed. The results of this experiment are shown in the histogram in Figure 14. The 
removal of 0743-6719 or 0302-2223 results in the two values of Y that are well-dehned and 
that are in reasonable agreement with the value calculated at high redshift in the constant 
threshold case. Removing only the one line from 0743-6719 nearest the Ly -a emission line 
with £ a bs = 1.5058 and observed equivalent width equal to 0.23 A results in Y = 6.23 x 10” 12 
s -1 . This object was part of the HST Key Project sample (Jannuzi et al. 1998) and they 
cite no evidence of associated aborption in its spectrum. Removing only the one line from 
0302-2223 nearest the Ly-a emission line with ^ a b s = 1.3886 and observed equivalent width 
equal to 0.27 A results in Y = 8.14 x 10 -12 s -1 . This object shows an absorption system 
at z a bs = 1.406 and is classified as an associated absorber. No metal absorption is seen 
at z a bs = 1-3886, though this absorber is within 5000 km s -1 of the QSO, the canonical 
associated absorber region. Removing both of these lines gives Y = 3.88 x 10~ 12 s -1 . Due 
to the small equivalent widths of both of these lines they are not included in the constant 
threshold analysis, and the solutions for J(v o) and Y for z > 1 are well-defined. 

It appears that this method has some trouble reliably recovering the background from a 
sample of absorption lines above an equivalent width threshold allowed to vary with S/N. As 
the method works well for the constant threshold case, we contend that the photoionization 
model, expressed in Equ. 3, used to create the likelihood function must not be an adequate 
model for the proximity effect when weak lines are included in the analysis. Liske & Williger 
(2001) introduce a method for extracting J{v 0 ) from QSO spectra based on flux statistics. 
We shall return to this topic in future work. 
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7. Discussion 
7.1. Radio Loudness 

As the results listed in Table 4 indicate, the inclusion of the four blazars and one BL 
Lac object, all at z < 1, in our sample does not change the result significantly. However, 
there is much observational evidence that radio loud and radio quiet quasars inhabit different 
environments, namely that radio loud quasars reside in rich clusters while radio quiet quasars 
exist in galaxy environments consistent with the held (Stockton 1982, Yee & Green 1984, 
1987, Yee 1987, Yates, Miller, & Peacock 1989, Ellingson, Yee, & Green 1991, Yee & Ellingson 
1993, Wold et al. 2000, Smith, Boyle, & Maddox 2000). If there is a corresponding increase 
in the number of Ly-cc absorption lines in the spectra of radio loud objects, this could cause 
the proximity effect to be suppressed, and the measured log[J(z/ 0 )] to be artificially large. 
We have therefore divided our sample into radio loud and radio quiet subsamples using the 
ratio of radio to UV flux to characterize the radio loudness, 

RL = log[S(5 GHz)]/log[S(1450 A)]. (15) 

The value of RL for each object in our sample is listed in Table 3. A histogram of these 
values and the distribution of RL with z for the sample objects are shown in Figure 15. The 
division between radio loud and radio quiet was chosen to be RL=1.0. The resulting values 
of log[J(z/ 0 )] for these subsamples are listed in Table 4. There is no significant trend for 
log[J(z/ 0 )] to appear larger for radio loud objects than for radio quiet objects. 


7.2. Non-Zero Ha 


We performed the maximum likelihood calculation for the case of a non-zero cosmolog¬ 
ical constant. This means that the observer-QSO and absorber-QSO luminosity distances 
that appear in the relationship between u and z (BDO) must be calculated numerically from 
the expression: 


d-L = 


(1 + % 


r dz' 

Jo W~Y 


(16) 


where 

E(z) = + Hk(l + z) 2 + Ha, (17) 


(Peebles, 1993) as this integral cannot be reduced to an analytical form for Ha A 0- 


The calculations in the sections above assume (H m ,Ha) = (1.0,0.0). Here, we perform 
the maximum likelihood search for J[y 0 ) using (Hm,Ha) = (0.3,0.7). For a QSO at z = 0.5 
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with a Lyman limit flux density of 0.1 /jJy, an absorber at z = 0.48, and an assumed 
background of log[J(z/ 0 )] = —22., this (O m ,^a) results in a value of u that is ~ 25% smaller 
than that inferred in the Ha = 0 case. Unlike all the other solutions performed, we ignore 
redshift path associated with metal lines and use all redshifts between z ( f- and z\ This 
does not change the results significantly, but cuts down the computation time substantially. 
The results are listed in Table 4 and are plotted in Figure 11. For comparison, we also give 
the solutions for J{y 0 ) found using the standard parameters, (H m ,Ha) = (1.0,0.0), with this 
redshift path neglected. We find that (Hm,Ha) = (0.3,0.7), does not change the value of 
J{v o) derived significantly from the value found using (Hm,Ha) = (1.0,0.0). 

We performed a slightly modified re-analysis of the Scott et al. (2000b) sample of ob¬ 
jects at z ~ 2 and found little effect at high redshift as well. The solution found for 
(H m ,Ha) = (1.0,0.0) was log[J(z/ 0 )]= —21.09^on7, while for (H m ,Ha) = (0.3,0.7), we find 
log[J(z/ 0 )]= —21.25lou7 for these data. 


7.3. dAf/dz 

In the case of a size distribution of Ly-a absorbers that is constant in redshift, the 
evolution of the number of Ly-cc absorption lines per unit redshift is given by: 

dAf/dz = A/”o(l + ^)“[HAf(l + ^) 3 + (1 — H.m — Ha)(1 + z ) 2 + Ha] °' 5 , (18) 

(Sargent et al. 1980) where Ao equals the absorber cross section times the absorber comov¬ 
ing number density times the Hubble distance, Trr^cfocHQ 1 . A plot of dAf/dz versus z for 
non-evolving Ly-cc absorbers in (fi M ,fi a) = (1.0,0.0) and (0.3,0.7) cosmologies is shown in 
Figure 17. It is clear that non-evolving models are too shallow to fit points at z > 1.7, so the 
normalization is found from a fit to the FOS data. The FOS data at z < 1.7 are consistent 
with a non-evolving population for (H m ,Ha) = (1.0, 0.0). The data are less consistent with 
a non-evolving concordance model in which (Hat, Ha) = (0.3, 0.7), though not significantly 
so. 


The number density evolution of Ly-a absorbers over the redshift range z = 0 — 5 cannot 
be approximated with a single power law. There is a significant break in the slope of the 
line number density with respect to redshift, near z = 1.7 (Weymann et al. 1998, Paper 
IV) though Kim, Cristiani, & D’Odorico (2001) argue that the break occurs at z = 1.2. 
Dave et al. (1999) show from hydrodynamical simulations of the low redshift Ly-ct forest, 
that the evolution of the line density is sensitive mainly to the HI photoionization rate, but 
also to the evolution of structure (cf. their Figure 7). The flattening of dAf/dz observed by 
Weymann et al. (1998) is mostly attributed to a dramatic decline in r(z) with decreasing z. 
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Dave et al. (1999) derive an expression for the density of Ly-a forest lines per unit redshift 
as a function of the HI photoionization rate: 

^ = C[(l + zfr-\z)f- 1 H-\z):, (19) 

where C is the normalization at some fiducial redshift which we choose to be z = 0 and T(z) 
can be expressed by Equ. 14. 

We fit the FOS and MMT absorption line data, binned in dAf/dz as presented in Paper 
IV and Scott et al. (2000a, Paper I), to this function in order to derive the parameters 
describing r(^) implied by the evolution in Ly-a forest line density. We observe flattening of 
dAf/dz at z < 1.7, but not to the degree seen by Weymann et al. (1998) in the Key Project 
data. As described in Paper IV, we find 7 = 0.54 ± 0 . 21 , for lines above a 0.24 A threshold, 
while Weymann et al. (1998) measure 7 = 0.15 ± 0.23. See Paper IV for more discussion 
of the significance and underlying causes of this difference. We find (Ahm, L>hm, z c , S ) = 
(3.0xl0“ 12 , 0.61, 5.5xl0 -7 , 7.07) and (lAxHT 11 , 0.38, 3.4xl0" 7 , 6.21) for (fi M , ^a) = (1-, 0.) 
and lines with rest equivalent widths above 0.24 and 0.32 A respectively. These fits to Equ. 19 
are shown in Figure 18(a). In panel (b), we plot r(z), as expressed in Equ. 14, evaluated 
using the parameters found from the fit to Equ. 19 above. The HM96 solution and the 
solution derived from the full FOS and MMT data sets are represented by the thick and 
thin solid lines respectively. The small values of z c derived from dAf/dz above translate 
into ionization rates that do not decrease dramatically with decreasing redshift and result 
from the less pronounced flattening of dAf/dz relative to the Key Project. These fits are 
particularly insensitive to the normalization, Ahm, so the errors on this parameter are large. 
These fits should therefore not be interpreted as measurements of r(z) as reliable as those 
found directly from the absorption line data. But we find them instructive nonetheless. The 
observed r(z) falls short of the ionization rate needed to fully account for the change in the 
Ly-a line density with redshift, indicating that if the value of 7 at low redshift is indeed 
slightly larger than that found by the Key Project, dAf/dz may still be consistent with a 
non-evolving population of Ly-a absorbers in the sense noted above, but the formation of 
structure in the low redshift universe must play a significant role in determining the character 
of the Ly-a forest line density. 


7.4. Comparison with Previous Results 

7.J h l. Proximity Effect 

KF93 performed a similar measurement with a small subsample of this total sample- 
the HST Quasar Absorption Line Key Project data of Bahcall et al. (1993). We compare 
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our result to that from Sample 2 of KF93, which was constructed from the Bahcall et al. 
(1993) data excluding one BAL quasar and all heavy element absorption systems. The Key 
Project sample has since been supplemented (Bahcall et al. 1996, Jannuzi et al. 1998) and 
those data have been included when appropriate in the complete archival sample of FOS 
spectra presented in Paper III. 

The mean intensity KF93 derive from their Sample 2 (6 = 35 km s” 1 , (3— 1.48, 7=0.21) 
is 5.0l 20 4 x 10” 24 ergs s” 1 cm -2 Hz” 1 sr” 1 . This result is lower than ours for z < 1 by a 
factor of ~ 13, though the errors are large on both results are large enough that they are 
consistent. We use 162 lines in our low redshift solution for J[y o), 65 more than KF93. 


7.4.2. Direct Measurements 

Several authors have examined the sharp cutoffs observed in the HI disks of galaxies 
in the context of using these signatures to infer the local ionizing background (Maloney 
1993, Corbelli & Salpeter 1993, Dove & Shull 1994). The truncations are modeled as arising 
primarily from photoionization of the disk gas by the local extragalactic background radiation 
held. Using 21 cm observations (Corbelli, Scheider, & Salpeter 1989, van Gorkom 1993) to 
constrain these models, limits on the local ionizing background are placed at 10 4 < <hi on < 
5 x 10 4 cm ” 2 s” 1 , where 

<F ion = 2n f fidfi [ ^-du = (20) 

J0 J u 0 ' W ' la s 

and where J v = I u for an isotropic radiation held. 

Additionally, narrow-band and Fabry-Perot observations of Ha emission from inter- 
galactic clouds (Stocke et al. 1991, Bland-Hawthorn et al. 1994, Vogel et al. 1995, Donahue, 
Aldering, & Stocke 1995) place limits of $ ion < 10 4 cm ” 2 s” 1 , or J(u 0 ) < 7.6 x 10 ” 23 ergs 
s ” 1 cm ” 2 Hz ” 1 sr ” 1 for a s = 1.8, while results from measurements of Galactic high velocity 
clouds (Kutyrev & Reynolds 1989, Songaila, Bryant, & Cowie 1989, Tufte, Reynolds, & 
Haffner 1998) imply <hi on < 6 x 10 4 cm ” 2 s” 1 , though the ionization of high velocity clouds 
may be contaminated by a Galactic stellar contribution. 

Tumlinson et al. (1999) have reanalyzed the 3C273/NGC3067 held using the Ha imaging 
data from Stocke et al. (1991) as well as new GHRS spectra of 3C273, in order to model 
the ionization balance in the absorbing gas in the halo of NGC3067. From this analysis, 
they derive the limits, 2600 < <bj 0 n < 10 4 cm ” 2 s” 1 , or 10 ” 23 < J(v 0 ) < 3.8 x 10 ” 23 ergs 
s ” 1 cm ” 2 Hz ” 1 sr ” 1 at z = 0.0047. Weymann et al. (2001) have recently reported an upper 
limit of < f > ion < 1.01 x 10 4 cm ” 2 s” 1 , or J{v 0 ) < 3.84 x 10 ” 23 ergs s ” 1 cm ” 2 Hz ” 1 sr ” 1 from 
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Fabry-Perot observations of the intergalactic HI cloud, 1225+01, for a face-on disk geometry. 
If an inclined disk geometry is assumed, this lower limit becomes J(z/ 0 ) < 9.6 x 1CU 24 ergs 
s^ 1 cm -2 Hz -1 sr” 1 . These results are summarized in Figure 16. It is encouraging that 
the proximity effect value is consistent with the limits on the background set by these more 
direct estimates which are possible locally. 


7.5. Comparison with Models 

Haardt & Madau (1996) calculated the spectrum of the UV background as a function of 
frequency and redshift using a model based on the integrated emission from QSOs alone. The 
QSO luminosity function is drawn from Pei (1995). The opacity of the intergalactic medium 
is computed from the observed redshift and column density distributions of Ly -a absorbers 
given by Equ. 2. The effects of attenuation and reemission of radiation by hydrogen and 
helium in Ly-a absorbers are included in these models. Their result for q 0 = 0.5 and a s = 1.8 
at z = 0 is J{v 0 ) = 1.6 x 10 -23 ergs s ” 1 cm -2 Hz ^ 1 sr _1 . 

Fardal et al. (1998) compute opacity models for the intergalactic medium (IGM) based 
on high resolution observations of the high redshift Ly-a: forest from several authors. Shull 
et al. (1999) extend the models of Fardal et al. (1998) to z = 0, treating opacity of low 
redshift Ly-a forest from observations made with HST/GHRS (Penton et al. 2000a,b) and 
with HST/FOS (Weymann et al. 1998). Like Haardt & Madau (1996), they also incorporate 
the observed redshift distribution of Lyman limit systems with log(Nni) > 17 (Stengler- 
Larrea et al. 1995, Storrie-Lombardi et al. 1994). Their models also allow for a contribution 
from star formation in galaxies in addition to AGN. The QSO luminosity function again 
is taken to follow the form given by Pei (1995) with upper/lower cutoffs at 0.01/10 L*. 
QSO UV spectral inclicies are assumed to equal 0.86, while the ionizing spectrum at v > u 0 
has a s = 1.8. The contribution to the background from stars was normalized to the Ha 
luminosity function observed by Gallego et al. (1995) and the escape fraction of photons of 
all energies from galaxies was taken to be < / esc >= 0.05. The full radiative transfer model 
described in Fardal et al. (1998) was used to calculate the contribution to the mean intensity 
by AGN, but not the contribution from stars, as they were assumed to contribute no flux 
above 4 Ryd, the energies at which the effects of IGM reprocessing become important. These 
authors find J{v 0 ) = 2.4 x lO -23 ergs s _1 cm -2 Hz -1 sr -1 at z ~ 0, with approximately equal 
contributions from AGN and stars, a value somewhat lower than our result for z < 1, but 
which is allowed within the errors. 

We estimate the contribution to the UV background from star-forming galaxies using 
the galaxy luminosity function of the Canada-France Redshift Survey (Lilly et al. 1995). At 
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z ~ 0.5, we derive J gal (Vo) = 1.5 x 10~ 22 ergs s _1 cm -2 Hz^ 1 sr -1 , assuming < f esc >= 1. The 
HM96 models for the QSO contribution give jQ so (z/ 0 ) = 5.2 x 10 -23 ergs s _1 cm -2 Hz -1 sr^ 1 
at z ~ 0.5. These estimates, and the range of measured J{y 0 ) in this paper, ~ 5 — 16 x 10~ 23 
ergs s _1 cm -2 Hz^ 1 sr _1 imply an escape fraction of UV photons from galaxies between 4% 
and 70%. The J{v o) inferred from dAf/dz in Section 7.3 implies escape fractions well over 
100 %. 

Bianchi et al. (2001) make updated estimates of the mean intensity of the background 
with contributions from both QSOs and star-forming galaxies. Their models incorporate 
various values of the escape fraction of Lyman continuum photons from galaxies which are 
constant with redshift and wavelength. Our new results at z < 1.7 are most consistent 
with their models of the QSO contribution alone, though some contribution from galaxies, 
ie. a small f esc , is allowed within the uncertainties. At z ~ 3.5, recent results from Steidel, 
Pettini, & Adelberger (2001) on the Lyman-continuum radiation from high redshift galaxies 
suggest that these sources become a more important component of the UV background at 
high redshift. 


7.6. Systematics 

Drawing on lessons learned from our work on high redshift objects in Paper II, we 
have made corrections for quasar systemic redshifts before performing the proximity effect 
analysis, as discussed in §3. This correction, ~ 300 km s -1 , was made to QSO redshifts 
measured from Ly-a emission for objects for which no systemic redshift measurement was 
available. For the low redshifts considered in this paper, redshifts measured from [OIII], 
Mgll, or Balmer emission lines were deemed suitable as QSO systemic redshift measurements. 

We have removed known gravitational lenses from the sample. As discussed above, we 
perform the proximity effect analysis omitting and including spectra that show associated 
absorption and damped Ly-a absorption and determined that neither of these populations 
significantly biases our results. 

Because we are working with low redshift data where line densities are low, we expect 
that blending has not contributed as strong a systematic effect as in the high redshift sample 
of Paper II. The curve-of-growth effects discussed in Paper II may still be present, since many 
lines in the sample have equivalent widths which place them on the flat part of the curve of 
growth. 

However, the effects of clustering may be even more important at low redshift than at 
high redshift. Loeb & Eisenstcin (1995) showed how the fact that quasars reside in the dark 
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matter potentials of galaxies and small groups of galaxies can influence the proximity effect 
signature. The peculiar velocities of matter clustered in these potentials can result in Ly-a 
absorption at redshifts greater than the quasar emission redshift. We found that including 
associated absorbers in our sample did not significantly change our results. Recently, Pas- 
carclle et ah (2001) report evidence for a lower incidence of Ly-a absorption lines arising in 
the gaseous halos of galaxies in the vicinities of QSOs than in regions far from QSOs. They 
argue that galaxy-QSO clustering may lead proximity effect measurements to overestimate 
J{y o) at z < 1 by a up to a factor of 20. While we agree that most systematic effects in this 
type of analysis, including clustering, will lead to overestimates of J{y 0 ), the agreement be¬ 
tween our results and the direct measurements discussed in Section 7.4.2 give us confidence 
that our results are not biased by this large a factor. 

The hydrodynamic simulations of the low redshift Ly-a forest of Dave et al. (1999) indi¬ 
cate that, at low redshift, structures of the same column density correspond to larger overden¬ 
sities and more advanced dynamical states than at high redshift. For a (Dm, Da) = (0.4, 0.6) 
cosmology, an equivalent width limit of 0.32 A corresponds to an overdensity of ~ 1.4 at 
z ~ 3, while at z ~ 0.6, this limit corresponds to Ph/J>h ~ 13. This may have implications 
on the clustering of Ly-a absorption lines around QSOs and hence on the values of J(z/ 0 ) 
derived from the proximity effect. It is possible that we are seeing this clustering effect in 
the variable threshold solution at z > 1 , in which the two highest c o(z) lines in the sample 
are responsible for the inability to isolate a reasonable maximum likelihood J{vq). 


8. Summary 

We have analyzed a set of 151 QSOs and 906 Ly-a absorption lines, the subset of the 
total data set presented in Paper III that is appropriate for the proximity effect. The primary 
results of this paper are as follows: 

(1) At low redshift, Balmer, [OIII], and Mg II emission lines are reasonable indicators 
of QSO systemic redshifts. Ly-a emission is blueshiftcd by ~ 300 km s ” 1 with respect to 

[OIII]. 

(2) The value of J{v 0 ) is observed to increase with redshift over the redshift range of the 
sample data, 0.03 < z < 1.67. Dividing the sample at z = 1, we find J{y 0 ) = 6.5l' 38 6 x 10” 23 
ergs s ” 1 cm -2 Hz ” 1 sr” 1 , at low redshift and J{v o) = I.OIq 8 x 10” 22 ergs s ” 1 cm ” 2 Hz ” 1 sr ” 1 
at high redshift. 

(3) The inclusion of blazars at z < 1 has no significant effect on the result. There is 
no significant difference between the values of J(// 0 ) derived from radio loud (RL > 1.0) and 
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radio quiet (RL < 1.0) objects, indicating that the observed richness of quasar environments 
does not distinctly bias the proximity effect analysis. 

(4) Using information measured and gathered from the literature on each QSO’s UV 
spectral index and solving for the HI ionization rate, yields 1.9 x 10 -13 s -1 for z < 1 
and 1.3 x 1CU 12 s _1 for and z > 1. Solving directly for the parameters (Ahm, Bum, Zc ’ &) 
in the HM96 parametrization of r(z) using the HST/FOS data presented by Bechtold et 
al. (2001) combined with the high redshift, ground-based data presented by Scott et al. 
(2000a,b) results in (Ahm, -Bhm, z c , S) = (7.6 x 1CT 13 , 0.35, 2.07, 1.77) for /3 = 1.46 and 
(A HM , H hm , z c , S ) = (3.2 x 10 " 13 , 1.45, 2.13, 1.42) for f3 = 1.7 for 1.7 < z < 3.8. 

(5) Allowing for a varying equivalent width threshold across each QSO spectrum results 
in consistently higher values of J{v 0 ) than are found from the constant threshold treatments. 
At z > 1, the variable threshold solution is not well-constrained. Jackknife experiments 
indicate that this is due the objects 0743-6719 and 0302-2223, namely the highest u(z) 
absorption lines in each of their spectra. 

( 6 ) Allowing for a cosmology in which (D M ,fbv) = (0.3, 0.7), rather than (l.,0.) has no 
significant effect on the value of J[y 0 ) derived from these data. 

(7) The z < 1 result is in agreement with the range of values of the mean intensity of the 
hydrogen-ionizing background allowed by a variety of local estimates, including Ha imaging 
and modeling of galaxy HI disk truncations. To within the uncertainty in the measurement, 
this result agrees with the one previous proximity effect measurement of the low redshift 
UV background (KF93). These results are consistent with calculated models based upon the 
integrated emission from QSOs alone (HM96) and with models which include both QSOs and 
starburst galaxies (Shull et al. 1999). The uncertainties do not make a distinction between 
these two models possible. 

( 8 ) The results presented here tentatively confirm the IGM evolution scenario provided 
by large scale hydrodynamic simulations (Dave et al. 1999). This scenario, which is success¬ 
ful in describing many observed properties of the low redshift IGM, is dependent upon an 
evolving J{v 0 ) which decreases from z = 2 to z = 0. However, the low redshift UV back¬ 
ground required to match the observations of the evolution of the Ly-a forest line density is 
larger than found from the data, indicating that structure formation is playing a role in this 
evolution as well. Our results and the work of others are summarized in Figure 16. We find 
some evidence of evolution in J[y 0 ), though it appears that even larger data sets, especially 
at z < 1 and/or improved proximity effect ionization models will be required to improve the 
significance. 
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Fig. 3.— Histograms of redshift differences between [OIII] and (a) Ly-a, (b) Mg II, and (c) 
Balmer emission lines, dotted lines show results from Laor et al. (1995) 
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Fig. 4.— FOS spectra used to extrapolate to the Lyman limit flux of each object: (dashed 
lines) power law fits, (dotted lines) la errors in spectrum 
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Remainder of Figure 4 available at: 

http: //lithops. as . arizona.edu/~jill/QuasarSpectra/ 

or 

http://hea-www.harvard.edu/QEDT/QuasarSpectra/ 
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Fig. 5.— Lyman limit luminosity versus redshift for objects in the HST/FOS sample 
(squares) and in the MMT sample presented in Papers I and II (crosses), solid line indi¬ 
cates the boundary between low and high luminosity objects 
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Fig. 6.— (a) Relative deficit of lines with respect to the number predicted by Equ. 1 for 
W t h. r = 0.32 A versus distance from the QSO for high and low luminosity QSOs (thick and 
thin solid lines, respectively) in both the HST/FOS sample presented in Paper III and the 
MMT sample presented in Paper I; (b) Deficit of lines within 2 h^ 1 Mpc as a function of QSO 
Lyman limit luminosity for the HST/FOS and MMT samples, the vertical line delineates 
the boundary between low and high luminosity objects 
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Fig. 7.— x 2 of binned data with respect to the ionization model expressed in Equ. 11 versus 
log[J(z/ 0 )] for various redshift ranges and equivalent width thresholds: (a) Wthr = 0.32 A; 
(b) W thr = 0.32 A, z <1; (c) W thr = 0.32 A, z > 1; (d) W thr = 0.24 A 
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Fig. 9.— Likelihood function versus log[J(z/ 0 )] for Wthr = 0.32 A (/3,b) = (a) (1.46,35); (b) 
(1.46,25); (c) (1.45,25); (d) (1.70,30); (e) (2.04,25); (f) (1.46,35), * < 1; (g) (1.46,35), z > 1; 
and for W thr = 0.24 A (/?, b) = (h) (1.46,35); (i) (1.46,25); (j) (1.45,25); (k) (1.70,30); (1) 
12.04.251 
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Fig. 10.— Number distribution per coevo 
J{y o) (KF method); (a-1) same as Fig. 9; 
divided by 5 for clarity 
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Fig. 11.— log[J(^o)] versus redshift, solid curves in (a)-(f) correspond to HM96 models: (a) 
W t hr = 0.32 A all redshifts, 0.03 < z < 1.67, and z < 1, z > 1 separately, ML method; (b) 
same as (a), BDO method; (c) variable threshold, all redshifts 0.03 < z < 1.67, and z < 1, 
z > 1 separately; (d) Wthr — 0-32 A all redshifts, RL > 0.3 and RL < 0.3; (e) Wthr — 0-32 A 
all redshifts, 0.03 < z < 1.67, and z < 1, total sample including blazars; (f) W t h r = 0.32 A 
z < 1 and z > 1, (solid points) (Q m ,^a) = (1.0, 0.0), (dotted points) — (0.3, 0.7), 

metal line dz neglected in both cases 
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Fig. 12.— (a) Values of log[J(z/o)] recovered from simulated QSO spectra with proximity 
effect included: (dotted lines)- input J{y 0 , z), see Figure 11(a) (solid points)- recovered J(u 0 ) 
for W t f ir = 0.32 A at all redshifts and at z < 1 and z > 1 separately; (b) same as (a), but 
J{v o) recovered using variable threshold 
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Fig. 13.— HI ionization rate versus redshift: (points)- constant equivalent width threshold 
maximum likelihood solutions from this paper, at z < 1 and z > 1, and from Paper II for 
1.7 < z < 3.8; (dashed line)- constant threshold solution to Equ. 13 for HST/FOS data 
alone; (solid line)- constant threshold solution to Equ. 14 with (3 = 1.46 and (3 = 1.7 for 
HST/FOS data and ground-based data from Papers I and II, (dotted line)- HM96 solution 
to Equ. 14 
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Fig. 14.— Histogram of results of jackknife measurements of HI ionization rate, T, for all 
lines at z > 1 above variable equivalent width threshold; labels on highest T bins indicate 
objects removed, see § 6.3 




Fig. 15.— (a) Histogram of radio loudness (RL) values for QSOs in proximity effect 
where RL = log[S(5 GHz)]/log[S(1450 A)], includes blazars and objects with damp 
absorption; (b) redshift versus RL for QSOs in proximity effect sample 





50 



log(1+z) 

Fig. 16.— log[ J(z/ 0 )] versus redshift: (lower limit at z ~ 0)- Tumlinson et al. (1999); (upper 
limit at z ~ 0)- Weymann et al. (2001); (filled triangle)- Shull et al. (1999); (upper limit 
at z — 0)- Weymann et al. (2001); (filled squares, bold error bars)- our results for z < 1 
and z > 1; (other filled squares)- results from KF93, Paper II, Lu et al. (1996), Savaglio et 
al. (1997), and Williger et al. (1994); (upper limit at z ~ 3)- Bunker et al. (1998); (solid 
curves)- HM96 models for two values of the global source spectral index, a s 





0 0.2 0.4 0.6 

log(1+z) 

Fig. 17.— (IN/dz versus z\ solid and dashed lines show the relation for non-evolving Ly-a 
absorbers given by Equ. 18 for (Ha^Ha) = (1.0, 0.0) and (0.3, 0.7), respectively; dotted lines 
are fits to low redshift data from Weymann et ah (1998) and to high redshift data of Kim 
et ah (1997); dashed-dotted lines are fits to low redshift data from Paper IV and to high 
redshift data from Paper I 
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1+z 



0 0.2 0.4 0.6 

log(1+z) 


Fig. 18.— (a) dAf/dz versus z: (solid points, dotted lines) Wthr = 0.24 A with fit to Equ. 19, 
(open points, dashed lines) W f y t/r = 0.32 A with fit to Equ. 19, (thick solid line) Equ. 19 
evaluated with HM96 parameters for r(z) expressed by Equ. 14, (thin solid lines) Equ. 19 
evaluated with parameters for r(z) found in this paper; (b) r(z) versus redshift expressed 
by Equ. 14 using HM96 parameters (thick solid line), using parameters found in this paper 
(thin solid lines), and using parameters found from fits to dAf/dz for W t j ir = 0.24 A and 
(Qm, Ov) = (1-0, 0.0) (dotted line) and W t hr = 0.32 A and Ry) = (1-0, 0.0) (dashed 
line), 








Table 1. Sample QSOs and 


QSO 1 NED description Ly-a Mgll OIII 

(a) (b) (c) 


0003+1553 opt.var. 
0003+1955 opt.var. 
0007+1041 opt.var. 
0015+1612 RQQ 
0017+0209 LINER 
0024+2225 • • • 

0026+1259 Syl 
0042+1010 
0043+0354 BAL? 2 
0044+0303 Syl? 
0050+1225 Compact,Syl 
0100+0205 opt.var. 
0102-2713 
0107-1537 
0112-0142 3 • • • 

0115+0242 3 opt.var. 
0117+2118 ••• 

0121-5903 Syl 
0122-0021 opt.var.,LPQ 
0137+0116 opt.var. 
0159-1147 opt.var.,Syl 
0214+1050 opt.var. 


0.4497 

0.4502 

0.4503 

0.0264 

0.0264 

0.0261 

0.0902 

0.0890 

0.089 

0.5492 



0.3994 



1.1081 

1.1096 


0.1453 

0.1463 

0.1452 

0.5854 

0.583 

0.586 

0.3803 



0.6219 

0.6222 


0.0594 



0.3937 


0.3936 

0.7763 



0.8574 



1.3739 

1.3727 


0.6652 

0.6700 


1.4925 

1.499 

1.504 

0.0461 

0.0462 

0.044 

1.0710 

1.0895 


0.2622 


0.2631 

0.6683 

0.6696 



0.4068 


0.407 


Line Redshifts 


Balmer _ References _ 

(d) (a) (b) (c) (d) 

••• (1) (2) (3) 

••• (1) (1) (4) 

0.0895 (1) (1) (5) (6) 

••• ( 1 ) 

••• ( 1 ) 

••• (1) (7) 

0.1458 (1) (1) (5) (6) 

0.584 (1) (8) (8) (8) 

••• ( 1 ) 

••• ( 1 ) ( 2 ) 

••• ( 1 ) 

••• (1) (3) 

••• ( 1 ) 

••• ( 1 ) 

••• ( 1 ) ( 1 ) 

••• (1) (9) 

1.499 (1) (10) (11) (11) 

••• (1) (1) (5) 

••• ( 1 ) ( 12 ) 

0.2644 (1) (1) (1) 

••• (1) (13) 

••• ( 1 ) ( 14 ) 



Table 1 



Balmer 


References 


(d) 

(a) 

(b) 

(c) 

(d) 


(1) 

(i) 




(1) 





(15) 





(1) 





(1) 

(7) 




(1) 




0.6206 

(1) 

(16) 


(a 


(1) 

(2) 




(1) 


(14) 


0.5731 

(1) 

(16) 

(14) 

(16) 


(1) 

(2) 

(5) 



(1) 

(13) 




(1) 





(1) 

(2) 

(14) 



(1) 

(10) 



0.0339 

(1) 

(1) 


(17) 


(1) 

(18) 



0.3698 

(1) 

(1) 

(1) 

(1) 

0.6570 

(1) 

(18) 


(18) 


(1) 




0.4620 

(1) 

(19) 

(14) 

(10) 

1.511 

(1) 

(20) 


(21) 





Table 1—Continued 


QSO 1 

NED description 

Ly-a 

(a) 

Mgll 

(b) 

OIII 

(c) 

Balmer 

(d) 

(a) 

References 
(b) (c) 

(d) 

0827+2421 

blazar,HPQ 

0.9363 

0.94 


0.942 

(1) 

(7) 


(7) 

0844+3456 

Syl 

0.0637 

0.0646 

0.064 


(1) 

(1) 

(5) 


0848+1623 

opt.var. 


1.9220 




(7) 



0850+4400 


0.5132 

0.5142 


0.5150 


(1) 


(1) 

0859-1403 3 

blazar 

1.3338 

1.3381 


1.341 

(1) 

(13) 


(21) 

0903+1658 3 

opt.var. 

0.4108 

0.4106 

0.4114 


(1) 

(22) 

(22) 


0907-0920 3 


0.630 4 








0916+5118 


0.5520 

0.5525 


0.5536 


(1) 


(1) 

0923+3915 3 

opt.var.,Syl,LPQ 

0.6986 

0.6990 



(1) 

(24) 



0935+4141 


1.937 4 








0945+4053 

LPQ 

1.2479 

1.2506 



(1) 

(19) 



0947+3940 

Syl 

0.2057 


0.2059 


(1) 


(25) 


0953+4129 

Syl? 

0.2331 


0.247 

0.2326 

(1) 


(25) 

(25) 

0954+5537 3 blazar,HPQ 

0.9005 

0.9025 



(1) 

(1) 



0955+3238 

opt.var.,Syl.8 

0.5281 


0.531 

0.5309 



(14) 

(10) 

0958+5509 


1.7569 

1.7582 



(10) 

(7) 



0959+6827 


0.7663 

0.7724 



(1) 

(1) 



1001+0527 


0.1589 

0.1605 


0.160 

(1) 

(1) 


(25) 

1001+2239 


0.9766 




(1) 




1001+2910 

AGN 

0.3285 



0.3293 

(1) 



(1) 

1007+4147 


0.6110 

0.6125 



(1) 

(13) 



1008+1319 


1.3012 

1.2968 



(1) 

(1) 








Table 1 


QSO 1 NED description Ly-ct Mgll OIII 

(a) (b) (c) 


1010+3606 Syl 
1026-004A • • • 

1026-004B 

1038+0625 opt.var.,LPQ 
1049-0035 Syl 
1055+2007 opt.var. 
1100+7715 opt.var.,AGN 
1104+1644 opt.var.,Syl 
1114+4429 Syl 
1115+4042 Syl 
1116+2135 E2,Syl? 
1118+1252 opt.var. 
1127-1432 3 blazar,LPQ 
1130+1108 ••• 

1136-1334 Syl 
1137+6604 opt.var.,LPQ 
1138+0204 • • • 

1148+5454 opt.var. 
1150+4947 opt.var. 
1156+2123 

1156+2931 blazar,HPQ 
1206+4557 • • • 


0.0785 


0.079 

1.4349 



1.5253 



1.2667 

1.272 


0.3580 

0.360 


1.1136 

1.1165 


0.3120 


0.324 

0.6294 


0.630 

0.1448 

0.1442 

0.143 

0.1545 

0.1552 




0.1768 

0.6823 



1.1824 

1.2121 


0.5065 


0.5110 

0.5551 

0.5571 


0.6449 

0.6448 

0.646 

0.3789 


0.3820 

0.9688 

0.9777 


0.3334 

0.333 

0.333 

0.3464 


0.3475 

0.7225 

0.7281 


1.1596 

1.164 



Balmer _References 


(d) 

(a) 

(b) 

(c) 

(d) 


(1) 


(5) 



(1) 





(1) 





(1) 

(7) 



0.3605 

(1) 

(5) 


(10) 


(1) 

(13) 



0.339 

(1) 


(25) 

(25) 

0.6307 

(1) 


(5) 

(6) 


(1) 

(1) 

(25) 


0.156 

(1) 

(1) 


(25) 

0.1756 



(25) 

(25) 


(1) 





(1) 

(18) 



0.5104 

(1) 


(1) 

(1) 

0.5604 

(1) 

(18) 


(18) 


(1) 

(13) 

(5) 


0.3831 

(1) 


(1) 

(1) 


(1) 

(10) 



0.333 

(1) 

(26) 

(26) 

(26) 

0.3459 

(1) 


(1) 

(1) 


(1) 

(1) 




(1) 

(7) 







Table 1 



Balmer 


References 


(d) 

(a) 

(b) 

(c) 

(d) 

0.0810 

(1) 

(i) 

(25) 

(25) 

0.3726 

(1) 



(1) 


(1) 




0.3374 

(1) 

(25) 

(5) 

(25) 


(1) 

(1) 




(1) 

(1) 

(5) 


0.158 

(1) 


(27) 

(27) 


(1) 

(13) 



0.4153 

(1) 


(1) 

(1) 


(1) 

(7) 




(10) 





(1) 

(28) 




(1) 

(7) 




(1) 





(1) 





(1) 




0.536 

(1) 

(29) 

(29) 

(29) 


(1) 

(7) 




(1) 




0.4853 

(1) 

(25) 


(25) 

0.2868 

(1) 

(12) 

(5) 

(6) 



Table 1—Continued 


QSO 1 

NED description 

Ly-a 

(a) 

Mgll 

(b) 

OIII 

(c) 

B aimer 

(d) 

(a) 

References 
(b) (c) 

(d) 

1305+0658 


0.6009 

0.5999 



(1) 

(1) 



1309+3531 

Sab,Syl 

0.1841 


0.184 

0.183 

(1) 


(25) 

(25) 

1317+2743 


1.0082 

1.016 



(1) 

(7) 



1317+5203 3 

blazar 

1.0550 

1.0555 



(1) 

(7) 



1318+2903 

opt.var. 

0.5469 




(1) 




1320+2925 


0.9601 

0.972 



(1) 

(7) 



1322+6557 

Syl 

0.1676 



0.1684 

(1) 



(25) 

1323+6530 


1.6227 

1.6233 



(1) 

(30) 



1327-2040 


1.1682 

1.170 



(1) 

(18) 



1328+3045 

DLAs 

0.8466 

0.8508 



(1) 

(13) 



1329+4117 


1.9351 




(10) 




1333+1740 


0.5464 

0.5546 



(1) 

(25) 



1351+3153 


1.3170 

1.3382 



(1) 

(31) 



1351+6400 

Syl 

0.0886 

0.0884 

0.087 

0.089 

(1) 

(1) 

(25) 

(25) 

1352+0106 


1.1200 




(1) 




1352+1819 

Syl 

0.1508 

0.1514 

0.1572 

0.1538 

(1) 

(1) 

(25) 

(25) 

1354+1933 

opt.var. 

0.7190 

0.718 

0.719 


(1) 

(7) 

(5) 


1356+5806 3 


1.3741 

1.370 



(1) 

(7) 



1401+0952 3 


0.4363 




(1) 




1404+2238 

Sy 

0.0966 

0.0978 


0.098 

(1) 

(1) 


(25) 

1407+2632 


0.95 

0.946 


0.958 

(1) 

(32) 


(32) 

1415+4509 


0.1145 

0.1142 

0.1143 

0.1139 

(1) 

(1) 

(25) 

(25) 






Table 1—Continued 


QSO 1 

NED description Ly-a 
(a) 

Mgll 

(b) 

OIII 

(c) 

B aimer 

(d) 

(a) 

References 
(b) (c) 

(d) 

1416+0642 


1.4339 



1.442 

(1) 



(21) 

1424-1150 


0.8033 

0.8037 



(1) 

(18) 



1425+2645 

opt.var. 

0.3634 



0.3644 

(1) 



(10) 

1427+4800 

Syl 

0.2215 


0.2203 

0.2246 

(1) 


(25) 

(25) 

1435-0134 


1.3099 




(1) 




1440+3539 

compact 

0.0764 

0.0772 

0.0777 

0.0772 

(1) 

(1) 

(25) 

(25) 

1444+4047 

El? 

0.2659 


0.2672 

0.267 

(1) 


(3) 

(5) 

1512+3701 

Syl? 

0.3704 

0.3734 

0.371 

0.3715 

(1) 

(2) 

(5) 

(6) 

1517+2356 


1.9037 




(10) 




1517+2357 


1.834 4 








1521+1009 


1.3210 

1.332 



(1) 

(7) 



1538+4745 


0.7704 

0.7711 



(1) 

(7) 



1544+4855 


0.3985 



0.4010 

(1) 



(2) 

1555+3313 3 


0.9402 

0.9427 



(1) 

(31) 



1611+3420 3 

blazar,LPQ 

1.3968 

1.3997 



(1) 

(33) 



1618+1743 

opt.var. 

0.5549 

0.5560 

0.555 


(1) 

(14) 

(13) 


1622+2352 

opt.var. 

0.9258 

0.925 



(1) 

(7) 



1626+5529 

Syl 

0.1315 

0.1325 

0.132 

0.133 

(1) 

(1) 

(25) 

(25) 

1630+3744 


1.4712 

1.478 

1.474 

1.478 

(1) 

(10) 

(11) 

(27) 

1634+7037 


1.3338 

1.338 

1.336 

1.342 

(1) 

(10) 

(11) 

(27) 

1637+5726 3 

LPQ 

0.7499 

0.750 


0.751 

(1) 

(7) 


(5) 

1641+3954 3 

opt.var.,HPQ 5 

0.5946 

0.5954 

0.593 


(1) 

(14) 

(2) 








Table 1—Continued 


QSO 1 

NED description 

Ly-a 

(a) 

Mgll 

(b) 

OIII 

(c) 

B aimer 

(d) 

(a) 

References 
(b) (c) 

(d) 

1704+6048 

opt.var. 

0.3694 

0.3704 

0.372 


(1) 

(2) 

(5) 


1715+5331 


1.9371 

1.932 



(10) 

(7) 



1718+4807 


1.0809 

1.0828 



(1) 

(7) 



1803+7827 

BL Lac 

0.6840 


0.6797 


(1) 



(23) 

1821+6419 

Syl 

0.2957 


0.297 


(1) 


(5) 


1845+7943 

opt.var.,BLRG,Syl 

0.0567 

0.0548 



(1) 

(1) 



2112+0556 


0.4585 



0.460 

(1) 



(5) 

2128-1220 

opt.var.,LPQ,Syl 

0.4988 

0.5000 

0.499 

0.5028 

(1) 

(2) 

(14) 

(6) 

2135-1446 

El,opt.var., Syl 

0.2016 


0.200 

0.199 

(1) 


(14) 

(34) 

2141+1730 

opt.var.,LPQ,Syl 

0.2124 


0.211 


(1) 


(14) 


2145+0643 

opt.var.,LPQ 

0.9997 

1.000 



(1) 

(7) 



2155-3027 3 

opt.var.,BL Lac 

0.116 4 








2201+3131 3 

LPQ 

0.2953 

0.2981 

0.295 

0.2979 

(1) 

(16) 

(5) 

(16) 

2216-0350 3 

opt.var. ,LPQ 

0.8997 

0.900 



(1) 

(7) 



2223-0512 3 

opt.var.,HPQ,BL Lac 

1.4037 




(1) 




2230+1128 3 

blazar,HPQ 

1.0367 

1.0379 



(1) 

(13) 



2243-1222 

opt.var.,HPQ 

0.6257 

0.6297 



(1) 

(17) 



2251+1120 

opt.var. 


0.322 

0.326 

0.3255 


(34) 

(5) 

(10) 

2251+1552 

blazar,HPQ 

0.8557 




(1) 




2251-1750 

opt.var., Syl 

0.0651 

0.0637 

0.064 


(1) 

(1) 

(5) 


2300-6823 


0.5149 

0.511 

0.516 

0.512 

(1) 

(35) 

(35) 

(35) 

2340-0339 


0.8948 

0.893 



(1) 

(7) 







Table 1—Continued 


QSO 1 

NED description 

Ly-a 

(a) 

Mgll 

(b) 

OIII 

(c) 

Balmer 

(d) 

(a) 

References 
(b) (c) 

(d) 

2344+0914 

opt.var.,Syl 

0.6710 

0.6722 

0.673 

0.6731 

(1) 

(16) (5) 

(16) 

2352-3414 

opt.var. 

0.7060 

0.7063 



(1) 

(2) 



1 See Paper III, Table 1 for alternate names 

2 We classify this as an associated absorber, see Paper III 

3 0bserved only with pre-COSTAR FOS and A-l aperture 

4 Redshift from Green et al. 1986 (0935+4141), Hewitt & Burbidge 1987 (1517+2357), Knezek & 
Bregman 1998 (0907-0920), Falomo et al. 1993 (2155-3027) 

5 Classified as blazar by Kinney et al. 1991 

References. — (l)This paper; (2)Tytler et al. 1987; (3)Stockton & MacKenty 1987; (4)de Robertis 
1985; (5)Corbin & Boroson 1996; (6)Zheng & Sulentic 1990; (7)Steidel & Sargent 1991; (8)Smith et al. 
1977; (9)Cristiani & Koehler 1987; (10)Tytlcr & Fan 1992; (ll)Nishihara et al. 1997; (12)Browne et 
al. 1975; (13)Aldcroft et al. 1994; (14)Corbin 1997; (15)Bolton et al. 1976; (16)Gaskcll 1982; (17)Basu 
1994; (18)Wilkes 1986; (19)Wills & Wills 1976; (20)di Serego-Alighieri et al. 1994; (21)Cheng et al. 
1990; (22)Lynds et al. 1966; (23)Lawrence et al. 1996; (24)Burbidge & Kinrnan 1966; (25)Green et 
al. 1986; (26)Lynds & Wills 1968; (27)Morris & Ward 1988; (28)Zotov 1985; (29)Netzer et al. 1994; 
(30)Barthel et al. 1990; (31)Ulrich 1976; (32)M c Dowell et al. 1995; (33)Schmidt 1977; (34)Kinman & 
Burbidge 1967; (35)Jauncey et al. 1978 



Table 2. Summary of Observations 


Name 

V 

Setup a 

Date 

Total exposure time (sec) 

0112-0142 

18.0 

1 

13Decl996 

1200 

0137+0116 

17.1 

1 

13Decl996 

1200 

0232-0415 

16.4 

1 

13Decl996 

1200 

0349-1438 

16.2 

1 

12Decl996 

900 

0414-0601 

15.9 

1 

19Decl995 

400 

0454-2203 

16.1 

1 

19Decl995 

400 

0624+6907 

14.2 

1 

19Decl995 

465 

0827+2421 

17.2 

3 

15Febl997 

1200 

0850+4400 

16.4 

1 

19Decl995 

300 

0859-1403 

16.6 

2a 

12Decl996 

3600 

0916+5118 

16.5 

1 

19Decl995 

350 

0923+3915 

17.9 

2b 

14Janl996 

1800 

0954+5537 

17.7 

2c 

20Aprl996 

3600 

0959+6827 

16.4 

2b 

14Janl996 

1800 

1001+2910 

15.5 

2a 

12Decl996 

3600 

1008+1319 

16.2 

2a 

10Decl996 

1800 

1130+1108 

16.9 

2d 

14Janl996 

3600 

1138+0204 

17.6 

2e 

12Decl996 

2400 

1156+2123 

17.5 

2e 

12Decl996 

1800 

1156+2931 

17.0 

2a 

10Decl996 

1800 

1214+1804 

17.5 

2f 

2lAprl996 

1800 

1230+0947 

16.1 

2f 

2lAprl996 

3600 

1305+0658 

17.0 

2c 

20Aprl996 

3600 


a Tclescope/Instrument set up: (1) FLWO 1.5 m, FAST, 300 1 mm 1 
1 st order, 3 " slit, 3660-7540 A; (2) SO B90, B&C, 600 1 mm -1 1 st order, 
1.5 " slit, [a] 4500-6700 A, [b] 3600-5825 A, [c] 4140-6370 A, [d] 6870- 
9140 A, [e] 5610-7860 A, [f] 5280-7550 A; (3) MMT, Blue Channel, 
800 1 mm - 1 1 st order, 2" slit, 4365-6665 A 







Table 3. Spectrophotometric Properties 


QSO 

N HI 
(a) 

robs 

(b) 

a 

(c) 

b 0 

(d) 

fobs 

(e) 

RL 

(f) 

(b) 

References 

(c) (d) 

(e) 

0003+1553 

3.88 

0.46 

0.71+0.52 

1.39+0.33 

1.94 (1450) 

2.24 

(2) 

(lb) 

(lb) 

(lb) 

0003+1955 

3.99 

2.04 

0.47+0.09 

6.77+0.45 

8.43 (1450) 

-0.44 

(3) 

(la) 

(la) 

(la) 

0007+1041 

5.62 


-0.50+1.00 

1.86+0.66 

1.47 (1450) 

0.00 


(la) 

(la) 

(la) 

0015+1612 

4.07 


-1.14+0.43 

0.19+0.06 

0.11 (1450) 

0.00 


(lb) 

(lb) 

(lb) 

0017+0209 

3.05 


1.98+0.56 

0.12+0.08 

0.31 (1450) 

0.00 


(lb) 

(lb) 

(lb) 

0024+2225 1 

3.60 


0.59+0.65 

0.60+0.15 

0.79 (1450) 

2.40 


(lc) 

(lc) 

(lc) 

0026+1259 

4.56 


-0.10+0.24 

2.33+0.33 

2.22 (1450) 

-0.04 


(lb) 

(lb) 

(lb) 

0042+1010 

5.52 


0.19+0.08 

0.08+0.02 

0.09 (1450) 

2.99 


(lc) 

(lc) 

(lc) 

0043+0354 1 

3.18 


2.35+0.04 

0.13+0.01 

0.97 (2093) 

0.00 


(lc) 

(lc) 

(lc) 

0044+0303 

2.88 

1.16 

0.34+0.11 

0.67+0.07 

0.79 (1450) 

1.94 

(2) 

(lc) 

(lc) 

(lc) 

0050+1225 1 

1.46 


0.84+1.14 

1.72+1.05 

2.56 (1450) 

0.06 


(la) 

(la) 

(la) 

0100+0205 

2.92 


1.42+0.27 

0.23+0.06 

0.45 (M50) 

0.00 


(lb) 

(lb) 

(lb) 

0102-2713 

1.93 



0.18 

0.29 (1285) 

0.00 


(4) 


(lb) 

0107-1537 

1.73 


0.78+0.31 

0.11+0.01 

0.16 (1450) 

0.00 


(lc) 

(lc) 

(lc) 

0112-0142 2 

4.32 



0.17 

0.29 (1326) 

3.83 


(4) 


(lc) 

0115+0242 2 

3.32 


0.83+0.08 

0.05+0.01 

0.08 (1450) 

4.08 


(lc) 

(lc) 

(lc) 

0117+2118 1 

4.75 

0.39 

0.15+40.6 

1.77+7.84 

1.88 (1307) 

0.00 

(2) 

(lc) 

(lc) 

(lc) 

0121-5903 

3.05 


0.15+0.10 

2.71+0.27 

2.91 (1450) 

0.00 


(la) 

(la) 

(la) 

0122-0021 

3.57 


0.65+0.07 

0.63+0.07 

0.86 (1450) 

3.13 


(lc) 

(lc) 

(lc) 

0137+0116 

3.00 


1.44+0.31 

0.03+0.02 

0.07 (1450) 

3.97 


(lb) 

(lb) 

(lb) 

0159-1147 

1.77 


-0.02+0.11 

1.35+0.05 

1.33 (1450) 

3.01 


(lc) 

(lc) 

(lc) 

0214+1050 

6.96 


1.39+0.08 

0.64+0.13 

1.22 (1450) 

2.57 


(lb) 

(lb) 

(lb) 

0232-0415 

2.42 

0.59 




2.73 

(2) 




0253-0138 2 

5.61 


0.31+0.19 

0.67+0.07 

0.78 (1450) 

0.00 


(lc) 

(lc) 

(lc) 

0254-3327B 4 

2.32 



0.28 


3.08 


(4) 



0302-2223 1,4 

1.87 

0.31 

-2.89+0.08 

2.57+0.44 

0.88 (1318) 

0.00 

(2) 

(lc) 

(lc) 

(lc) 

0333+3208 

13.5 


0.80+5.79 

0.56+1.65 

0.81 (1450) 

3.38 


(lc) 

(lc) 

(lc) 

0334-3617 2 

1.40 


0.13+1.27 

0.13+0.02 

0.14 (1450) 



(lc) 

(lc) 

(lc) 

0349-1438 

3.83 


-0.32+0.29 

2.45+0.23 

2.11 (1450) 

2.53 


(lc) 

(lc) 

(lc) 

0355-4820 1 

1.16 

0.39 

0.65+0.58 

0.52+0.13 

0.70 (1450) 

2.91 

(5) 

(lc) 

(lc) 

(lc) 






Table 3—Continued 


QSO 

N hi 
(a) 

robs 

i ^o 

(b) 

a 

(c) 

b 0 

(d) 

fobs 

(e) 

RL 

(f) 

(b) 

References 

(c) (d) 

(e) 

0403-1316 2 

3.65 


0.23+0.04 

0.35+0.05 

0.39 (1450) 

4.35 


(lc) 

(lc) 

(lc) 

0405-1219 

3.74 

2.05 

-0.11+0.04 

4.38+0.18 

4.14 (1450) 

2.68 

(2) 

(lc) 

(lc) 

(lc) 

0414-0601 

5.14 

0.34 

-0.19+0.08 

0.77+0.05 

0.70 (1450) 

2.66 

(2) 

(lc) 

(lc) 

(lc) 

0420-0127 3 

7.10 


1.84+0.05 

0.08+0.01 

0.20 (1450) 

3.89 


(lc) 

(lc) 

(lc) 

0439-4319 

2.30 


0.40+0.08 

0.27+0.01 

0.33 (1450) 

2.95 


(lc) 

(lc) 

(lc) 

0454+0356 5 

7.39 

0.38 

-0.26+2.26 

1.40+0.57 

1.26 (1336) 

2.50 

(2) 

(lc) 

(lc) 

(lc) 

0454-2203 

2.99 

0.38 

0.05+4.19 

1.25+0.17 

1.28 (1450) 

2.77 


(lb) 

(lb) 

(lb) 

0518-4549 

4.12 


0.18+1.45 

0.12+0.05 

0.13 (1450) 

5.06 


(la) 

(la) 

(la) 

0537-4406 2 

4.02 

0.05 

2.00+0.16 

0.14+0.03 

0.36 (1450) 

4.05 

(2) 

(lc) 

(lc) 

(lc) 

0624+6907 

7.01 


1.71+0.03 

2.37+0.18 

5.26 (1450) 

0.00 


(lb) 

(lb) 

(lb) 

0637-7513 

9.22 

0.53 

1.32+0.07 

0.27+0.03 

0.49 (1450) 

4.10 

(2) 

(lc) 

(lc) 

(lc) 

0710+1151 2 

11.0 


0.16+0.08 

1.13+0.10 

1.22 (1450) 

4.12 


(lc) 

(lc) 

(lc) 

0742+3150 1 

4.89 

0.35 

0.24+0.43 

0.92+0.08 

1.03 (1450) 

2.96 

(2) 

(lb) 

(lb) 

(lb) 

0743-6719 

11.9 

0.24 




3.46 

(2) 




0827+2421 3 

3.51 


1.21+0.04 

0.34+0.03 

0.59 (1450) 

3.17 


(lc) 

(lc) 

(lc) 

0844+3456 1 

3.31 


0.75+0.03 

2.31+0.09 

4.94 (2495) 

0.00 


(lc) 

(lc) 

(lc) 

0848+1623 4 

29.7 


0.46 

0.15 

0.19 (1450) 

0.00 


(6) 


(11) 

0850+4400 

2.53 


1.02+0.20 

0.35+0.05 

0.56 (1450) 

0.00 


(lb) 

(lb) 

(lb) 

0859-1403 2 

5.71 

0.60 




3.29 

(2) 




0903+1658 2 

3.61 


3.28+0.27 

0.03+0.02 

0.17 (1450) 

2.79 


(lb) 

(lb) 

(lb) 

0907-0920 6 

4.57 


-0.04+1.50 

0.11+0.008 

0.11 (1822) 

0.00 


(lc) 

(lc) 

(lc) 

0916+5118 

1.40 


0.31+0.03 

0.71+0.06 

0.82 (1450) 

0.00 


(lc) 

(lc) 

(lc) 

0923+3915 2 

1.53 


0.17+0.05 

0.70+0.03 

0.77 (1450) 

4.83 


(lc) 

(lc) 

(lc) 

0935+4141 4,E 

’1.32 



0.55 


0.00 


(4) 



0945+4053 

1.44 


-0.33+5.03 

0.17+0.19 

0.15 (1450) 

4.07 


(lc) 

(lc) 

(lc) 

0947+3940 

1.61 


0.70+0.11 

0.90+0.08 

1.25 (1450) 

0.00 


(lb) 

(lb) 

(lb) 

0953+4129 

1.28 


0.71+0.08 

1.13+0.10 

1.58 (1450) 

0.10 


(lb) 

(lb) 

(lb) 

0954+5537 2 

0.94 


0.96+0.05 

0.12+0.01 

0.18 (1450) 

3.51 


(lc) 

(lc) 

(lc) 

0955+3238 1 

1.62 

0.38 

0.96+0.07 

0.45+0.03 

0.87 (1774) 

2.99 

(2) 

(lc) 

(lc) 

(lc) 

0958+5509 1 

0.84 

0.31 




0.00 

(2) 









Table 3—Continued 


QSO 

N hi 

robs 

L vo 

(b) 

a 

b 0 

fobs 

^1/ 

RL 


References 


(a) 

(c) 

(d) 

(e) 

(f) 

(b) 

(c) 

(d) 

(e) 

0959+6827 

3.93 


1.12+2.21 

0.54+0.71 

1.10 (1720) 

1.99 


(lc) 

(lc) 

(lc) 

1001+0527 1 

2.41 


1.73+0.12 

0.24+0.04 

0.55 (1450) 

0.26 


(lb) 

(lb) 

(lb) 

1001+2239 

2.82 


1.67+0.32 

0.05+0.02 

0.12 (1450) 

3.17 


(lc) 

(lc) 

(lc) 

1001+2910 

1.93 


1.18+0.02 

1.08+0.06 

1.88 (1450) 

0.00 


(lb) 

(lb) 

(lb) 

1007+4147 

1.23 

0.72 

-0.20+0.08 

1.12+0.07 

1.02 (1450) 

2.92 

(2) 

(lc) 

(lc) 

(lc) 

1008+1319 1 

3.79 

0.58 




0.00 

(2) 




1010+3606 

1.24 


0.90+1.60 

0.66+0.60 

1.00 (1450) 

0.00 


(la) 

(la) 

(la) 

1026-004A 

4.85 



0.11 

0.19 (1328) 

0.00 


(4) 


(lc) 

1026-004B 

4.85 



0.15 

0.24 (1285) 

0.00 


(4) 


(lc) 

1038+0625 1 

2.81 


-0.65+1.96 

1.30+0.06 

1.00 (1361) 

3.09 


(lc) 

(lc) 

(lc) 

1049-0035 1 

3.87 

0.35 

1.60+0.11 

0.51+0.07 

1.07 (1450) 

0.00 

(2) 

(lb) 

(lb) 

(lb) 

1055+2007 

1.94 


0.51+0.37 

0.27+0.05 

0.34 (1450) 

3.64 


(lc) 

(lc) 

(lc) 

1100+7715 1 

3.04 


0.67+0.04 

0.97+0.06 

1.33 (1450) 

2.76 


(lb) 

(lb) 

(lb) 

1104+1644 

1.55 


-0.02+0.15 

1.23+0.08 

1.22 (1450) 

2.66 


(lc) 

(lc) 

(lc) 

1114+4429 1 

1.80 


1.80+0.04 

0.15+0.02 

0.35 (1450) 

0.00 


(lb) 

(lb) 

(lb) 

1115+4042 1 

1.86 


0.44+0.05 

1.10+0.14 

1.35 (1450) 

0.00 


(lb) 

(lb) 

(lb) 

1116+2135 

1.27 


0.46+0.10 

2.31+0.36 

2.87 (1450) 

0.01 


(lb) 

(lb) 

(lb) 

1118+1252 1 

2.28 


0.42+0.08 

0.11+0.02 

0.14 (1450) 

2.75 


(lc) 

(lc) 

(lc) 

1127-1432 2 

4.07 


0.96+3.07 

0.31+0.59 

0.49 (1450) 

4.78 


(lc) 

(lc) 

(lc) 

1130+1108 1 

3.47 


1.40+0.15 

0.32+0.05 

0.62 (1450) 

0.00 


(lb) 

(lb) 

(lb) 

1136-1334 

3.51 

0.60 

-0.46+0.20 

1.03+0.09 

0.83 (1450) 

3.36 

(2) 

(lb) 

(lb) 

(lb) 

1137+6604 1 

1.00 

1.05 

0.24+0.04 

1.04+0.09 

1.17 (1450) 

2.98 

(2) 

(lc) 

(lc) 

(lc) 

1138+0204 1 

2.37 


0.97+0.09 

0.22+0.03 

0.35 (1450) 

0.00 


(lb) 

(lb) 

(lb) 

1148+5454 

1.19 

0.97 

0.56+0.17 

1.04+0.11 

1.35 (1450) 

-0.13 

(2) 

(lc) 

(lc) 

(lc) 

1150+4947 

2.01 


0.66+0.05 

0.19+0.03 

0.26 (1450) 

3.44 


(lb) 

(lb) 

(lb) 

1156+2123 

2.56 


0.95+0.10 

0.31+0.04 

0.49 (1450) 

2.23 


(lb) 

(lb) 

(lb) 

1156+2931 

1.58 

0.57 

1.27+0.08 

0.73+0.06 

1.33 (1450) 

3.04 

(2) 

(lc) 

(lc) 

(lc) 

1206+4557 

1.27 

0.45 

-0.32+0.49 

1.96+0.23 

1.69 (1450) 

0.00 

(2) 

(lc) 

(lc) 

(lc) 

1211+1419 

2.70 


1.27+0.34 

1.31+0.37 

2.37 (1450) 

-0.37 


(la) 

(la) 

(la) 

1214+1804 1 

2.74 


1.55+0.17 

0.25+0.05 

0.52 (1450) 

0.00 


(lb) 

(lb) 

(lb) 






Table 3—Continued 


QSO 

N hi 
(a) 

robs 

i ^o 

(b) 

a 

(c) 

b 0 

(d) 

fobs 

(e) 

RL 

(f) 

(b) 

References 

(c) (d) 

(e) 

1215+6423 1 

2.10 


-0.14+2.50 

0.19+0.06 

0.18 (1340) 

3.18 


(lc) 

(lc) 

(lc) 

1216+0655 

1.57 


0.84+0.06 

0.97+0.06 

1.44 (1450) 

0.44 


(lb) 

(lb) 

(lb) 

1216+503a 7 

1.87 



0.35 

0.58 (1326) 

0.00 


(4) 

(lc) 

(lc) 

1219+0447 1 

1.68 


0.83+0.05 

0.06+0.006 

0.15 (2457) 

0.00 


(lc) 

(lc) 

(lc) 

1219+7535 2 

3.13 


0.00+0.36 

2.21+0.34 

2.21 (1450) 

0.45 


(la) 

(la) 

(la) 

1226+0219 3 

1.81 

7.40 

-1.51+2.68 

47.6+1.94 

26.9 (1330) 

4.26 

(7) 

(la) 

(la) 

(la) 

1229-0207 5 

2.34 

0.23 

1.23+0.78 

0.32+0.20 

0.57 (1450) 

3.25 


(lc) 

(lc) 

(lc) 

1230+0947 2 

1.81 


1.33+0.36 

0.51+0.16 

0.96 (1450) 

0.00 


(lb) 

(lb) 

(lb) 

1241+1737 

1.81 

0.25 




2.16 

(2) 




1247+2647 5 

1.03 

0.76 




-0.07 

(2) 




1248+3032 

1.23 


0.19+0.28 

0.08+0.01 

0.09 (1450) 

3.19 


(lc) 

(lc) 

(lc) 

1248+3142 8 

1.27 



0.26 


0.00 


(4) 


(8) 

1248+4007 

1.44 

0.57 

0.67+0.76 

0.48+0.16 

0.65 (1450) 

0.00 

(2) 

(lc) 

(lc) 

(lc) 

1249+2929 8 

1.14 



0.22 


0.00 


(4) 


(8) 

1250+3122 

1.24 



0.33 

0.54 (1279) 

0.00 


(4) 

(lb) 

(lb) 

1252+1157 

2.34 


0.80+0.38 

0.37+0.07 

0.54 (1450) 

3.12 


(lc) 

(lc) 

(lc) 

1253-0531 2 

2.12 

1.43 

1.58+0.02 

0.14+0.01 

0.30 (1450) 

4.47 

(2) 

(lc) 

(lc) 

(lc) 

1257+3439 1 

1.13 



0.51 

0.94 (1450) 

1.14 


(4) 


(9) 

1258+2835 1 

0.93 


0.21+0.81 

0.32+0.04 

0.34 (1331) 

0.00 


(lc) 

(lc) 

(lc) 

1259+5918 

1.37 

1.02 

1.14+0.45 

0.96+0.32 

1.63 (1450) 

0.00 

(2) 

(lb) 

(lb) 

(lb) 

1302-1017 

3.37 

0.99 

1.17+0.06 

2.00+0.14 

3.47 (1450) 

2.34 

(2) 

(lb) 

(lb) 

(lb) 

1305+0658 

2.16 


-0.07+0.04 

0.24+0.03 

0.23 (1450) 

3.13 


(lc) 

(lc) 

(lc) 

1309+3531 1 

2.55 


1.08+0.16 

0.68+0.11 

1.12 (1450) 

1.58 


(lb) 

(lb) 

(lb) 

1317+2743 

1.18 

0.73 

0.64+0.19 

1.04+0.10 

1.40 (1450) 

0.00 

(2) 

(lc) 

(lc) 

(lc) 

1317+5203 1,2 

! 1.90 


0.54+0.82 

0.51+0.15 

0.66 (1450) 

2.70 


(lc) 

(lc) 

(lc) 

1318+2903 

1.14 

0.26 

-0.06+10.0 

0.58+0.25 

0.56 (1450) 

0.00 

(2) 

(lb) 

(lb) 

(lb) 

1320+2925 

1.17 


1.37+1.63 

0.19+0.26 

0.36 (1450) 

0.00 


(lc) 

(lc) 

(lc) 

1322+6557 

1.92 


0.91+0.16 

0.66+0.07 

1.01 (1450) 

0.00 


(lb) 

(lb) 

(lb) 

1323+6530 1,4, f .99 



0.11 


3.02 


(4) 



1327-2040 1 

7.53 

0.19 

0.83+0.41 

0.55+0.12 

0.82 (1450) 

2.62 

(2) 

(lc) 

(lc) 

(lc) 






Table 3—Continued 


QSO 

N hi 
(a) 

robs 

L vo 

(b) 

a 

(c) 

b 0 

(d) 

fobs 

(e) 

RL 

(f) 

(b) 

References 

(c) (d) 

(e) 

1328+3045 5 

1.16 


0.39+0.13 

0.20+0.01 

0.24 (1450) 

4.49 


(lc) 

(lc) 

(lc) 

1329+4117 5 

0.97 

0.95 




0.00 

(2) 




1333+1740 

1.75 

0.51 

0.92+4.71 

0.65+1.81 

1.01 (1450) 

1.39 

(2) 

(lb) 

(lb) 

(lb) 

1351+3153 1 - 

’1.29 


-0.91+2.78 

0.16+0.28 

0.11 (1319) 

2.88 


(lc) 

(lc) 

(lc) 

1351+6400 1 

2.10 


0.97+0.06 

1.62+0.07 

4.36 (2531) 

1.10 


(lc) 

(lc) 

(lc) 

1352+0106 

2.25 

0.07 

0.50+1.23 

0.83+0.33 

1.05 (1450) 

0.00 

(2) 

(lc) 

(lc) 

(lc) 

1352+1819 

2.03 


0.38+0.13 

0.59+0.11 

0.71 (1450) 

0.00 


(lb) 

(lb) 

(lb) 

1354+1933 

2.21 

0.40 

0.68+0.11 

0.56+0.05 

0.77 (1450) 

3.53 

(2) 

(lc) 

(lc) 

(lc) 

1356+5806 2 

1.40 


0.09+6.29 

0.56+0.04 

0.59 (1344) 

2.34 


(lc) 

(lc) 

(lc) 

1401+0952 2 

1.96 


2.03+0.29 

0.12+0.05 

0.31 (1450) 

0.72 


(lb) 

(lb) 

(lb) 

1404+2238 1 

1.99 


1.03+0.04 

0.31+0.04 

0.86 (2413) 

0.29 


(lc) 

(lc) 

(lc) 

1407+2632 

1.47 

0.83 

0.28+0.05 

1.20+0.07 

1.38 (1450) 

0.00 

(2) 

(lc) 

(lc) 

(lc) 

1415+4509 

1.13 


0.65+0.08 

0.85+0.05 

1.32 (1790) 

0.00 


(lb) 

(lb) 

(lb) 

1416+0642 1 

6.24 


1.20+15.9 

0.25+3.47 

0.40 (1308) 

3.67 


(lc) 

(lc) 

(lc) 

1424-1150 

7.54 


-0.04+0.18 

0.85+0.07 

0.83 (1450) 

2.59 


(lc) 

(lc) 

(lc) 

1425+2645 1 

2.55 

0.15 

1.67+0.10 

0.22+0.03 

0.48 (1450) 

2.43 

(2) 

(lb) 

(lb) 

(lb) 

1427+4800 

1.88 


0.47+0.24 

0.69+0.07 

0.86 (1450) 

0.00 


(lb) 

(lb) 

(lb) 

1435-0134 

3.66 

0.82 




0.00 

(5) 




1440+3539 1 

1.00 


0.44+0.09 

3.61+0.15 

4.96 (1857) 

-0.58 


(lb) 

(lb) 

(lb) 

1444+4047 

1.27 

0.89 

0.86+0.04 

1.06+0.06 

1.59 (1450) 

0.00 

(2) 

(lb) 

(lb) 

(lb) 

1512+3701 

1.39 

0.57 

0.94+0.12 

0.61+0.06 

0.95 (1450) 

2.75 

(2) 

(lb) 

(lb) 

(lb) 

1517+2356 4 

3.91 



0.51 


0.00 


(4) 



1517+2357 4 

3.91 



0.08 


0.00 


(4) 



1521+1009 

2.88 

1.65 




0.00 

(2) 




1538+4745 1 

1.64 

0.34 

0.57+0.06 

1.03+0.05 

1.34 (1450) 

1.28 

(2) 

(lc) 

(lc) 

(lc) 

1544+4855 

1.60 

0.10 

2.04+1.72 

0.36+0.81 

0.95 (1450) 

0.00 

(2) 

(lb) 

(lb) 

(lb) 

1555+3313 2 

2.35 


1.79+0.08 

0.03+0.005 

0.07 (1450) 

3.03 


(lc) 

(lc) 

(lc) 

1611+3420 2 

1.65 



0.18 

0.30 (1322) 

4.88 


(4) 

(lc) 

(lc) 

1618+1743 

4.14 


-0.30+0.05 

1.30+0.06 

1.13 (1450) 

2.70 


(lb) 

(lb) 

(lb) 

1622+2352 

4.46 


1.75+0.16 

0.09+0.01 

0.21 (1450) 

3.54 


(lc) 

(lc) 

(lc) 






Table 3—Continued 


QSO 

(a) 

fobs 

+0 

(b) 

a 

(c) 

f,o 

(d) 

fobs 

(e) 

RL 

(f) 

(b) 

References 

(c) (d) 

(e) 

1626+5529 1 

1.83 


0.29+0.16 

1.13+0.22 

1.30 (1450) 

0.00 


(lb) 

(lb) 

(lb) 

1630+3744 

1.07 

0.84 




0.00 

(2) 




1634+7037 

4.55 

1.96 




0.00 

(2) 




1637+5726 2 

1.90 


0.17+0.02 

0.64+0.05 

0.70 (1450) 

3.98 


(lc) 

(lc) 

(lc) 

1641+3954 3 

1.02 

0.61 

1.04+0.08 

0.41+0.06 

0.67 (1450) 

3.92 

(2) 

(lc) 

(lc) 

(lc) 

1704+6048 1 

2.32 

0.90 

1.25+0.16 

0.94+0.14 

1.68 (1450) 

2.86 

(2) 

(lb) 

(lb) 

(lb) 

1715+5331 

2.71 


0.43 

0.58 

0.29 (1450) 

0.53 


(10) 


(2) 

1718+4807 1 

2.27 


-0.43+0.84 

5.01+1.16 

4.09 (1450) 

1.52 


(lc) 

(lc) 

(lc) 

1803+7827 3 

3.92 


1.69+0.02 

0.53+0.05 

1.16 (1450) 

3.35 


(lc) 

(lc) 

(lc) 

1821+6419 1 

3.98 

1.86 

0.86+0.07 

3.90+0.13 

8.37 (2204) 

1.10 

(2) 

(lc) 

(lc) 

(lc) 

1845+7943 1 

4.17 


0.66+0.27 

0.42+0.08 

0.58 (1450) 

3.88 


(la) 

(la) 

(la) 

2112+0556 

6.48 

0.29 

0.48+0.93 

0.54+0.16 

0.67 (1450) 

0.00 

(2) 

(lb) 

(lb) 

(lb) 

2128-1220 

4.75 

0.35 

0.27+1.63 

1.77+0.51 

2.02 (1450) 

2.99 

(11) 

(lb) 

(lb) 

(lb) 

2135-1446 1 

4.71 


0.94+0.42 

0.57+0.15 

0.88 (1450) 

3.17 


(lb) 

(lb) 

(lb) 

2141+1730 1 

8.20 


1.22+0.05 

0.81+0.14 

1.43 (1450) 

2.84 


(lb) 

(lb) 

(lb) 

2145+0643 

4.90 


0.99+0.68 

0.72+0.31 

1.14 (1450) 

3.58 


(lc) 

(lc) 

(lc) 

2201+3131 2 

9.02 

0.60 

0.96+0.08 

3.15+0.33 

4.93 (1450) 

3.64 

(2) 

(lb) 

(lb) 

(lb) 

2216-0350 2 

5.66 

0.18 

1.21+0.09 

0.40+0.05 

0.71 (1450) 

3.43 

(2) 

(lc) 

(lc) 

(lc) 

2223-0512 2 

5.47 

0.16 




4.35 

(2) 




2230+1128 2 

5.42 


0.76+0.96 

0.45+0.21 

0.64 (1450) 

4.39 


(lc) 

(lc) 

(lc) 

2243-1222 

4.94 


-0.38+0.06 

1.50+0.10 

1.25 (1450) 

3.32 


(lc) 

(lc) 

(lc) 

2251+1120 1 

5.08 


1.2 

1.46 

0.49 (1450) 

3.06 


(12) 


(2) 

2251+1552 

6.38 

0.09 

1.04+0.05 

0.71+0.07 

1.15 (1450) 

3.94 

(2) 

(lc) 

(lc) 

(lc) 

2251-1750 1 

2.77 


1.06+0.08 

1.47+0.09 

4.32 (2507) 

0.07 


(lc) 

(lc) 

(lc) 

2300-6823 

3.69 


-0.34+0.75 

0.26+0.04 

0.22 (1450) 

3.18 


(lb) 

(lb) 

(lb) 

2340-0339 

3.61 


0.68+0.05 

0.99+0.06 

1.36 (1450) 

2.24 


(lc) 

(lc) 

(lc) 

2344+0914 

5.76 

0.34 


0.22 

0.41 (1450) 

3.52 

(2) 

(4) 

(9) 

(9) 

2352-3414 

1.08 


0.07+0.03 

0.74+0.05 

0.77 (1450) 

2.70 


(lc) 

(lc) 

(lc) 


Note. — Columns (a) 10 20 cm -2 from Stark et al. 1992; Burstein & Heiles 1982; Lockman & 
Dickey 1995; (b) Direct measurement of flux at Lyman Limit in ynJy; (c) Measured spectral index; 
(d) Extrapolated flux at Lyman Limit in yuJy; (e) Measured flux at the rest wavelength listed in 
A; (f) Radio Loudness, RL = log[S(5 GHz)]/log[S(1450 A)] 

References. — (1) this paper, FOS data [a] H130, [b] H190, [c] H270; (2)Lanzetta et al. 
1993; (3)Zheng et al. 1995; (4)Zheng et al. 1997; (5)Hamann et al. 1995; (6) Tytler & Fan 1992; 
171 Anncnzcller et al. 1998: l8lSa,ndulea,k Pesch 1984: l9iOsrner et al. 1994: ilfll Zhcne - Malkan 







Table 4. Measurements of J[v 0 ) 


Sample 

Nlines 

7 ,norm. 

P 

b 

method 

log [</+o)] 

x 2 

Qx 2 

Qks 

(a) 

(b) 

(c) 


(d) 

(e) 

(f) 

(g) 

(h) 

(i) 

1 

259 

0.82, 13.6 

1.46 


BDO 

-22.04++ 

2.13 

0.95 

0.80 

1 

259 

0.82, 6.73 

1.46 

35 

ML 

-22 11 +0 ' 51 
ZZ.ll_o 40 

1.21 

0.29 

0.80 

1 

259 

0.82, 9.61 

1.46 

25 

ML 

no 4 O+0.52 

-ZZ.lZ_o.39 

1.01 

0.41 

0.80 

1 

259 

0.82, 9.31 

1.45 

25 

ML 

99 i q+0.51 

-ZZ.lD_o .44 

0.78 

0.58 

0.80 

1 

259 

0.82, 11.8 

1.70 

30 

ML 

_21 74 + 0 - 45 
Zi - ‘^-0.36 

1.34 

0.23 

0.80 

1 

259 

0.82, 38.0 

2.04 

25 

ML 

-21 47 + 0 - 43 
-0.32 

1.10 

0.35 

0.80 

2 

289 

0.15, 31.3 

1.46 


BDO 

-22 06 +ao5 
zz.uu_o 62 

2.62 

0.91 

0.30 

2 

289 

0.15, 12.0 

1.46 

35 

ML 

-22.03++ 

1.32 

0.24 

0.30 

2 

289 

0.15, 13.9 

1.46 

25 

ML 

-22.04+°;|l 

1.34 

0.23 

0.30 

2 

289 

0.15, 13.6 

1.45 

25 

ML 

-22.06+;$ 

1.48 

0.18 

0.30 

2 

289 

0.15, 17.6 

1.70 

30 

ML 

-21-69++ 

1.47 

0.18 

0.30 

2 

289 

0.15, 31.1 

2.04 

25 

ML 

-21 42 +0 - 37 

zi.+z_ 0 28 

0.88 

0.50 

0.30 

la 

162 

1.50, 10.1 

1.46 


BDO 

99 07+1-19 
-zz.cw _ a82 

1.51 

0.98 

0.64 

la 

162 

1.50, 4.92 

1.46 

35 

ML 

- 22 . 181 °;“ 

0.17 

0.98 

0.64 

la 

162 

1.50, 3.67 

1.46 

35 

ML 

91 79+1-52 1 
-zi. ( z_ 0 74 

1.02 

0.40 

0.62 

la 

162 

1.50, 3.71 

1.46 

35 

ML 

9 1 00+1-54 2 
-zl.OO _ 0 73 

0.98 

0.43 

0.62 

lb 

97 

-0.87, 53.0 

1.46 


BDO 

- 22 . 02++ 5 

2.44 

0.87 

0.98 

lb 

97 

-0.87, 26.1 

1.46 

35 

ML 

-21.98+54 

2.25 

0.03 

0.98 

lb 

97 

-0.87, 21.5 

1.46 

35 

ML 

-21.7618;“ 1 

1.31 

0.24 

0.95 

lb 

97 

-0.87, 21.5 

1.46 

35 

ML 

21 Q'S 4 ' 0 - 93 2 
zi.yo_Q.57 

1.27 

0.26 

0.95 

3 

214 

0.28, 9.97 

1.46 

35 

ML 

91 P- 7 +O. 8 O 
-Zl.O 1 _o . 5 2 

0.47 

0.82 

0.70 

4 

208 

1.04, 5.76 

1.46 

35 

ML 

99 1 c;+0.66 
-ZZ.iO_o.46 

1.47 

0.19 

0.65 

5 

373 

0.60, 7.93 

1.46 

35 

ML 

_21 74 + 0 - 55 
- 61 - ^—0.39 

0.97 

0.44 

0.96 

6 

301 

0.89, 6.57 

1.46 

35 

ML 

99 1 7+0-44 
-zz.l f —0.37 

0.98 

0.43 

0.97 

7 

415 

0.67, 7.72 

1.46 

35 

ML 

91 OO+0-46 
-Zl.OZ_ 0.37 

0.93 

0.46 

0.98 

7a 

213 

0.79, 7.28 

1.46 

35 

ML 

99 99 +O .74 
-ZZ.ZZ_o .54 

0.29 

0.94 

0.64 

7b 

202 

0.72, 7.29 

1.46 

35 

ML 

-21.60+;!? 

1.15 

0.33 

0.98 

8 

422 

0.69, 7.64 

1.46 

35 

ML 

-21 85 +0 ' 46 

Zl.00_o.36 

0.82 

0.55 

0.97 

8 a 

220 

0.84, 7.10 

1.46 

35 

ML 

99 oq+0.73 
-ZZ.ZO_o. 4 9 

0.46 

0.83 

0.56 

8 b 

202 

0.72, 7.29 

1.46 

35 

ML 

-21.60+;!? 

1.15 

0.33 

0.98 
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Table 4—Continued 


Sample 

Nlines 

7 ,norm. 

P 

b 

method 

log[J(z/ 0 )] 

x 1 2 

Q x 2 

Qks 

(a) 

(b) 

(c) 


(d) 

(e) 

(f) 

(g) 

(h) 

(i) 

9 

906 

0.61, 9.26 

1.46 

35 

ML 

91 91 +0.49 
-Z1.Z1 _ 0 32 

0.55 

0.76 

0.91 

9a 

474 

0.63, 9.23 

1.46 

35 

ML 

-21 7Q +0 - 53 

' a -0.40 

0.76 

0.59 

0.87 

9b 

432 

1.05, 6.40 

1.46 

35 

ML 

-20.821“, 

0.33 

0.91 

0.71 


Note. — 

(a) Sample number- ( 1 ) all lines with W > 0.32 A, (la) z < 1 , (lb) z > 1 ; ( 2 ) all lines 
with W > 0.24 A, (2a) z < 1, (2b) z > 1; (3) RL > 1.0; (4) RL < 1.0; (5) sample (1) 
including associated absorbers; ( 6 ) sample (1) including damped Ly-a absorbers; (7) sample 
(1) including both associated absorbers and damped Ly-a absorbers; (7a) z < 1, (7b) z > 1; 
( 8 ) sample (1) including associated absorbers, damped Ly-a absorbers, and blazars; ( 8 a) 
z < 1, ( 8 b) z > 1; (9) sample ( 8 ), all lines above variable threshold, (9a) z < 1, (9b) z > 1 

(b) Number of Ly-a forest lines in sample 

(c) Normalization, definition depends on method listed in (e), see § 6 , Equ. 9 and Equ. 10 

(d) Doppler parameter in km s ^ 1 

(e) BDO- Bajtlik, Duncan, & Ostriker (1988), ML- maximum likelihood, see Kulkarni & Fall 
(1993) 

(f) Best fit value of log[J+o)] hr units of ergs s _1 cm -2 Hz -1 sr -1 

(g) X 2 °f data versus the ionization model used 

(h) x 2 probability for the ionization model used 

(i) K-S probability 

1 no dz removed for metal lines 

2 no dz removed for metal lines, (9 m, ^a) = (0.3, 0.7) 







Table 5. Simulation Results 


Input log[J(z/ 0 )] 

2 w thr 

lA 

Recovered log[ J(z/ 0 )] X 2 

(a) 

(b) (c) 

(d) 

(e) (f) 


-23.0 


all 

0.32 

1.41,7.81 

-22.74«;i| 

5.13 

-23.0 


all 

variable 

1.15,8.94 

-22.477;!! 

5.73 

- 22.0 


all 

0.32 

1.17,8.27 

91 qq+0.41 
-Zl.OZ _ 0 29 

13.6 

- 22.0 


z < 1 

0.32 

0.95,8.74 

- 20 . 817 s 

9.53 

- 22.0 


z > 1 

0.32 

1.79,5.21 

-21.64«;! 

6.70 

- 22.0 


all 

variable 

1.48,6.71 

-21 63 +a20 
zi.oo _ 0 19 

3.21 

- 22.0 


z <1 

variable 

0.75,10.0 

-21.3418;“ 

11.3 

- 22.0 


z > 1 

variable 

1.52,6.09 

-21.631°;^ 

1.30 

- 21.0 


all 

0.32 

1.44,7.25 

-2o.8ii8;4i 

1.56 

- 21.0 


all 

variable 

1.13,8.72 

- 2 o. 8 ii 8 :l 

0.73 

(0.017) log(l + z) 

- 21.87 

all 

0.32 

0.99,9.46 

-21.5418:1 

3.35 

(0.017) log(l + z) 

- 21.87 

z < 1 

0.32 

0.51,11.1 

- 2 i. 8 oi 8 :l 

4.63 

(0.017) log(l T z') 

- 21.87 

z > 1 

0.32 

1.90,4.74 

-21.541818 

1.55 

(0.017) log(l + z) 

- 21.87 

all 

variable 

1.38,7.32 

-21.5618:1 

0.57 

(0.017) log(l + z) 

- 21.87 

z < 1 

variable 

0.84,10.1 

-21.8318:1 

3.28 

(0.017) log(l T z^j 

- 21.87 

z > 1 

variable 

2.48,2.67 

91 Q7+0-31 
-0.23 

0.82 


Note. — (a) Value of log[J(z/ 0 )] used for modifying absorber column densities according 
to Equ. 4 and Equ. 5 ; (b) Redshift range of solution; (c) Equivalent width threshold in A 
for line sample used in solution; (d) Maximum likelihood 7 for line sample used, maximum 
likelihood method normalization, see § 6 , Equ. 10; (e) Value of log[ J(z/ 0 )] from simulated 
spectra using the ML technique; (f) y 2 of data versus the BDO ionization model 







Table 6 . HI Ionization Rates 


Sample 7 ,A 

(a) (b) 

P 

b 

log [Thi] 

xl 

(c) 

Qxi 

(d) 

log[J(i/ 0 )] 

(e) 

1 

0.69,7.65 

1.46 

35 

19 1 7+0-50 
-iz.l /_q 40 

0.49 

0.81 

-21.56 

la 

0.85,7.11 

1.46 

35 

-12.7018;” 

0.38 

0.88 

-22.09 

lb 

0.72,7.29 

1.46 

35 

-ii.88i8;5o 

0.48 

0.81 

-21.28 

2 

0.61,9.27 

1.46 

35 

-11 27 +0 ' 74 
ii - z mo.45 

0.78 

0.58 

-20.67 

2 a 

0.63,9.24 

1.46 

35 

1 9 oq+0.55 
- 1 Z.ZO _0 42 

1.17 

0.31 

-21.62 

2 b 

1.05,6.40 

1.46 

35 

-9.089l 2 °° 22 

1.17 

0.31 

-18.48 

1 

0.69, 7.21 1 

1.46 

35 

-12.67,1.73 1 

1.01 

0.40 


2 

0.61, 9.04 1 

1.46 

35 

-10.86,3.04 1 

0.47 

0.82 



Note. — (a) ( 1 ) all lines with W > 0.32 A, (la) z < 1 , (lb) z > 1 , ( 2 )all 
lines above variable threshold, (2a) z < 1, (2b) z > 1; (b) Maximum likelihood 
method normalization (see § 6 , Equ. 10); (c) y 2 of data versus the ionization 
model used; (d) y 2 probability for the ionization model used; (e) J{y 0 ) implied 
by T listed and a s = 1.8 (see §6.2, Equ. 12) 

1 maximum likelihood solution for log(H p i),.Bpi and normalization (see § 6 . 2 , 
Equ. 13) 







